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Abstract 



We studied the phase structures of Af = 1 supersymmetric SO(N c ) gauge theory with 
Nf flavors in the vector representation as we deformed the M = 2 supersymmetric QCD by 
adding the superpotential of arbitrary polynomial for the adjoint chiral scalar field. Using weak 
and strong coupling analyses, we determined the most general factorization forms for various 
breaking patterns. We observed all kinds of smooth transitions for quartic superpotential. 



1 Introduction and summary 



The M = 1 supersymmetric gauge theories in four dimensions have rich structures and one 
obtains nonperturbative aspects by studying the holomorphic effective superpotential which 
determines the quantum moduli space. A large class of interesting gauge theories can be 
obtained from the choice of geometries in which D5-branes wrap partially over the nontrivial 
cycles [1, 2, 3, 4]. As far as the effective superpotential is concerned, the geometry in which 
the four dimensional gauge theories are realized on the worldvolume of D5-branes wrapping 
around S 2 is replaced by a dual geometry in which D5-branes are replaced by RR fluxes and 
the S 2 by S 3 . These RR fluxes provide the effective superpotential which corresponds to one 
of four dimensional gauge theories on the D5-branes [5, 6, 7]. This equivalence has been tested 
for several different models in different contexts [8, 9, 10, 11, 12, 13]. 

A new recipe for the computation of the effective superpotential was proposed by Dijk- 
graaf and Vafa [14, 15, 16] through the free energies in a certain matrix model. This matrix 
model analysis could be interpreted within purely field theoretical point of view without a 
string theory [17]. Based on the anomalous Ward identity of a generalized Konishi anomaly, 
a powerful machinery for the nonperturbative aspects was obtained. Moreover, a new kind of 
duality where one can transit several vacua with different broken gauge groups continuously 
and holomorphically by changing the parameters of the superpotential was found in [18]. The 
extension of [18] to the M = 1 supersymmetric gauge theories with the gauge group SO/USp 
was found in [19] in which there were no flavors and the phase structures of these theories, a 
matrix model curve and a generalized Konishi anomaly equation, were obtained. More recently 
[20], by adding the flavors in the fundamental representation to the theory of [18], the vacuum 
structures in classically and quantum mechanically were described and an addition map as well 
as multiplication map were developed. 

In this paper, we study the phase structures of M = 1 supersymmetric SO(N c ) gauge theory 
with Nf flavors in the vector representation by deforming the M = 2 supersymmetric QCD 
with the superpotential of arbitrary polynomial for adjoint chiral scalar field, by applying the 
methods in [18, 19, 20]. These kinds of study were initiated by Cachazo, Seiberg and Witten in 
[18], where a kind of new duality was found. This paper is a generalization of [19] to the case 
with flavors. We found that with flavors the phase structure is richer and that more interesting 
dualities show up. We list some partial papers [21]- [63] on the recent works, along the line of 
[14, 15, 16]. 

In section 2, we describe the classical moduli space of M = 2 SQCD deformed to M = 1 
theory by adding the superpotential W(<&) (2.2). The gauge group SO(N c ) will break to 
SO(N ) x n"=i U(Nj) with N + E"=i 2N j = N c by choosing the adjoint chiral field $ to be the 
root of W'(x) or ±mj, the mass parameters. To have pure Coulomb branch where no factor 
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U(Nj) is higgsed, we restrict ourselves to the case W'(±mi) = 0. For each factor with some 
effective massless flavors, there exists a rich structure of Higgs branches, characterized by an 
integer r i: that meets the Coulomb branch along the submanifold called the root of the r-th 
Higgs branch. 

In section 3.1, we discuss the quantum moduli space of SO(N) by both the weak and strong 
coupling analyses. When the difference between the roots of W'(x) is much larger than M = 2 
dynamical scale A (in other words, in the weak coupling region), the adjoint scalar field <3> 
can be integrated out and it gives a low energy effective J\f = 1 superpotential. Under this 
condition, one can suppress higher order terms except the quadratic piece in the superpotential 
(2.2). Then the effective superpotential consists of the classical part plus nonperturbative part. 
There exist two groups of solutions, i.e., Chebyshev vacua and Special vacua, according to the 
unbroken flavor symmetry, meson matrix M (3.1) and various phases of vacua. 

At the scales below the Af = 1 scale A x (when the roots of W'(x) are almost the same), the 
strong coupling analysis is relevant. We need to look for the special points where some number 
of magnetic monopoles (mutually local or non local) become massless, on the submanifold of the 
Coulomb branch of Af = 2 SO(N) which is not lifted by the AT — 1 deformation. The conditions 
for these special points are translated into a particular factorization form of the corresponding 
Seiberg-Witten curve. We discuss the characters of these curves at the Chebyshev branch or 
the Special branch, especially, the power of factor t = x 2 and the number of single roots. 

In section 3.2, combining the quantum moduli space of SO(N) group with the quantum 
moduli space of U (N) group studied in [20], we give the most general factorization forms of the 
curves with the proper number of single roots and double roots for various symmetry breaking 
patterns, which generalize the results in [9, 10, 11, 12]. From the point of view of the geometry, 
these various breaking patterns correspond to the various distributions of wrapping D5-branes 
among the roots of W'(x), i.e., some of the roots do not have wrapping D5-branes. We also 
summarize the counting of vacua for various phases which will be used heavily in the examples 
in sections 4 and 5. For U(Ni) without flavors, the number of vacua is given by the Witten 
index w(Ni) = JVj. For U(Ni) with flavors, it is given by (3.11) (for more details, see [20]). For 
SO(Ni) without flavors, it is also given by the Witten index and finally for SO(Ni) with 
flavors, the number of vacua is given by (iVj — Mj — 2) for the Chebyshev vacua and one for 
Special vacua. 

In section 4, we analyze the simplest nontrivial examples for quartic tree level superpotential 
with massive flavors for 5*0(4), 5*0(5), and 50(6) gauge groups 1 . For these cases, there are 
two values which can be chosen by $ and we have the following breaking patterns: SO(N c ) — > 

1 Although for the pure case [19], the description for SO(7) and SO (8) gauge theories was given, in this paper 
we will skip it because on the one hand, it will give rise to the complicated solutions and it is hard to analyze 
and on the other hand, we expect one cannot see any new interesting phenomena. 
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SO (N ) x UiNi), SO(N c ) -> SO(N c ), and 50(iV c ) -> Z7([iV c /2]). Depending on both the 
SO(N) factor at the Chebyshev branch or the Special branch and the U(N) factor at the various 
r baryonic or non-baryonic branches, the curve takes a different factorization form. By solving 
the factorization of the curve, we find various phases predicted by weak and strong coupling 
analyses that lead to the matches of the counting of vacua. We find also three interesting smooth 
transitions among these three breaking patterns: SO(N c ) <-> SO(M ) x [/(Mi), SO(N ) x 
C/(iVi) <-> 50 (M ) x C/(Mi), and £/([iV c /2]) <-> 50 (M ) x C/(Mi). Especially the smooth 
transitions of 50(iV c ) <-> 50 (M ) x [/(Mi) and Z7([JV c /2]) <-> 50 (M ) x [/(Mi) cannot be 
found if we do not consider the breaking patterns SO(N c ) -> [/([Ay2]) and SO(N c ) -> SO(N c ). 
We also discuss carefully how the smooth transition arises when we tune the parameters of the 
deformed superpotential. The phase structures for various gauge groups have been summarized 
in the Tables which can be found in sections 4 and 5. 

In section 5, we move to the massless flavors with quartic deformed superpotential for 
SO(N c ) where N c = 4,5,6 and 7. In this case, contrary to section 5, at the IR limit, the 
SO(N ) factor has massless flavors instead of a U(Ni) factor having effective massless flavors 
as in section 5. Because of this difference, new features arise. For example, there is Special 
branch for general SO(N c ) with massless flavors while for massive case, only 50(2), 50(3), 
and 50(4) gauge groups have the Special branch. Also with massless flavors, there are no 
smooth transitions as in SO(N c ) <-> SO(N ) x C/(7Vi) and 50(iV c ) <-> U([N c /2\). Furthermore, 
for the smooth transition 5O(iV ) x U(N 1 ) <->■ 5O(M ) x U(Mi) in the Special branch, we have 
M = (2Nf - No + 4), which is the relationship between SO(N ) and 50 (M ) to be the Seiberg 
dual pair. In fact, the smooth transition in this case may be rooted in the Seiberg duality. In 
this section, we also give the general discussion for the possible smooth transition in the case 
of massless flavors 2 . 

In Appendix A, by using the Af = 2 curve together with monopole constraints we are 
interested in and applying the contour integral formula, we derive the matrix model curve 
(A. 19) for deformed superpotential with an arbitrary degree and the relationship (A. 6) (or the 
most general expression (A. 8)) between the matrix model curve (i.e., the single-root part of the 
factorized Seiberg- Witten curve) and the deformed superpotential W'(x). The formula (A. 8) 
will be used to determine the number of vacua for fixed tree level superpotential. Using these 
results, we have also checked the generalized Konishi anomaly equation for our gauge theory 
with flavors (A. 24). 

In Appendix B, in order to understand the vacua of different theories, we discuss both the 
addition map and multiplication map. The addition map with massive flavors relates the vacua 
of SO(N c ) gauge theory with Nf flavors in the r-th branch to those of SO(N' c ) gauge theory 
2 In fact, with a little modification, the result can also be applied to the massive flavors. 
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with Njr flavors in the r'-branch where both N c and N' c are even or odd. This phenomenon is 
the same exactly as the one in the U(N C ) gauge theory with flavors. However, the addition 
map with massless flavors can relate SO(2N c ) gauge theory to SO(2N c + 1) gauge theory. 
For the multiplication map, we present the most general form. Through this multiplication 
map, one can obtain the unknown factorization of the gauge group with higher rank from the 
known factorization of the gauge group with lower rank. In this derivation, the properties 
of Chebyshev polynomials are used. In several specific examples, we demonstrate the general 
results. Some interesting points of these general multiplication maps are: (1) We can map the 
results of SO(2M) to SO(2N + 1) and vice versa; (2) We can map the results without massless 
flavors to the case with some massless flavors. This is a somewhat unusual result. 

One can also study the phase structures of M = 1 supersymmetric USp(2N c ) gauge theory 
with Nf flavors in the fundamental representation with equal footing described in this paper. 
However, due to the length of the paper, we will publish our results in a separate paper. 

Finally let us mention a few interesting directions to be done in the future: 

• In this paper, we restrict ourselves to the case W'(±m) = 0. It would be interesting to 
study the general case without this constraint W'(±m) = along the line of [50]. For the 
general case, the parameter space is bigger and we expect to have more smooth transitions in 
this bigger space. 

• We have discussed some indices which distinguish between different phases. It will be 
useful to have more concrete discussion of phases with the complete list of indices, along the 
line of [18]. 

• In this paper, we have discussed the phase structure with adjoint chiral field $. It would 
be interesting to discuss the phase structure with other representations, especially the second 
rank tensor representation. The descriptions with these representations have attracted some 
attention recently [47, 54, 55, 60]. 

2 The classical moduli space of SO(N c ) supersymmetric 
QCD 

In this section, we will discuss the classical moduli space of SO(N c ) with flavors and the 
Af = 1 deformation (2.2) so that the total superpotential will be (2.1). Although the classical 
picture will be modified by quantum corrections, it does give some useful information about 
the quantum moduli, especially in the weak coupling region. 

Let us consider an Af = 1 supersymmetric SO(N c ) gauge theory with Nj flavors of quark 
Q % a (i = 1, 2, • • • , 2Nf, a = 1, 2, • • • , N c ) in the vector representation. The tree level superpoten- 
tial of the theory is obtained from Af = 2 SQCD by adding the arbitrary polynomial of the 
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adjoint scalar $ ab belonging to the J\f = 2 vector multiplet: 

fc+i 

W tree ($, Q) = V^Qi^QiJij + V2m t] Q\Qi + E ^Tr$ 2s (2.1) 

where the first two terms come from the M = 2 theory and the third term, W(3>), can be 
described as a small perturbation of A/" = 2 SO(N c ) gauge theory [64, 65, 66, 67, 68, 69, 70, 
71, 72, 8, 11] 

fc+1 1 

WW = E f 1 ™ 25 = E fe« 2s , « 2s = -Tr$ 2s (2.2) 

s=l 2S s=l ZS 

where $ a b is an adjoint scalar chiral superfield that plays the role of a deformation breaking 
M = 2 supersymmetry to M = 1 supersymmetry. Note that in [64] the coefficient for the 
mass term of quark in the tree superpotential is different from a/2 used here, but it can be 
absorbed in the mass matrix and the only quadratic mass deformation for <3> with gi = fi (other 
parameters are vanishing) was considered in [64] . The is the symplectic metric used to raise 
and lower the SO(N c ) flavor indices while is a quark mass matrix, and they are defined as 

J = ^ ^ g ^ <E> f-N f xN f , m= ^l ) ® dia S( m l'"'' m ^/) ( 2 - 3 ) 

where Ijv/xjv> is the NfXNf identity matrix. The M = 2 theory without W^^) is asymptotically 
free for Nf < N c — 2 (and generates M = 2 strong-coupling scale A), conformal for Nf = N c — 2 
(scale invariant) and infrared (IR) free for Nf > N c — 2. 

The classical vacuum structure, the zeroes of the scalar potential, can be obtained by solving 
D-terms and F-terms. We summarize the following results from the mechanism of adjoint vevs 
as follows: 

1) The eigenvalues of $, ±i0j, can only be the roots of W'(x) or ±mj; thus, the gauge 
group SO(N c ) with Nf flavors is broken to the product of blocks with or without the effective 
massless flavors. Among these blocks at most one block is SO(/i n+ i) and others, U([/,i) where 
i — 1, 2, • • • , n. 

2) If (pi = rkirrii but W'(±mj) 7^ (W'(x) has some roots which are not equal to ±m,), 
the corresponding gauge symmetry of that block will be completely higgsed; Because of this, 
we will restrict the following discussions to the case of W'(±m) = for which there are richer 
structures. 

3) For each block, if it has "effective" massless flavors, there are various Higgs branches 
classified by some integers r. 
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3 The quantum moduli space of SO(N c ) supersymmetric 
QCD 

From the above section we see that the gauge group SO(N c ) is broken to the product 
of various U(Ni) with at most one SO(N). To get the quantum moduli space we need first 
to understand the quantum moduli for various factors. The quantum theory of U(N i ) with 
effective massless flavors has been discussed in [73, 74, 20] and we will not repeat it here. Here 
we will focus on the quantum moduli of SO(N) 3 with massless flavors discussed in [68, 64]. 

3.1 The quantum theory of SO(N) supersymmetric QCD 

The quantum theory of SO(N) with mass deformation |/zTr<3> 2 has been presented in [68, 64] 
for both weak and strong coupling analyses. Notice that when N < 4, it is not a simple Lie 
group, so we will restrict our discussion to the case N > 5 4 . 

3.1.1 The weak coupling analysis 

Special corners in parameter space will place M = 1 vacua in regions where the gauge 
symmetry breaking scale is much larger than the M = 2 dynamical scale, A. Then the gauge 
coupling is small. When the mass [i for adjoint scalar $ is larger compared with the dynamical 
scale A, the adjoint scalar $ can be integrated out to lead to a low energy effective M = 1 
superpotential from which one can recover the original theory as in [75]. The weak coupling 
analysis is valid when the difference between the roots of the polynomial W (x) are much bigger 
than A. Under this condition we can only consider the quadratic part of the effective superpo- 
tential since the higher powers of it will be suppressed by \i. Then the relevant superpotential 
is 

W tree ($,Q) = V2Q\$ ab Q{J l3 + ^Tr$ 2 . 

Integrating out $ will give 

W tree (<Z>, Q) = "^ Tr (M J M J) , M ij = Q i ■ Q j . (3.1) 

Then the effective superpotential will consist of the above classical part (3.1) plus nonpertur- 

bative effects. To find Af = I vacua, the effective superpotential should be minimized. We can 

3 The number of colors N is less than or equal to N c and it can be 2, 4, • • • , N c for N c even and 3, 5, • • • , N c 
for N c odd. 

4 For N > 5 the Witten index w(N) of SO(N) gauge group is (N — 2) while for N < 4 the Witten indices 
are 1, 1, 1, 2, 4 for N = 0, 1, 2, 3, 4 respectively. We will give some explanations why the Witten index is different 
for N < 4. 
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study the vacuum structure, the number of vacua, and global symmetry breakings according 
to the range of the number of flavors Nf. 

We summarize two groups of solutions as follows: 

1) The first group has non-degenerated meson matrix M. They are all in a confining phase 
with unbroken flavor symmetry U(Nf) and the number of vacua is given by (iV — Nf — 2). As we 
will show, they are given by Chebyshev point in the corresponding Seiberg-Witten(SW)-curve. 
Its position in the M = 2 moduli space is located by the roots of the Chebyshev polynomial; the 
light degrees of freedom are mutually nonlocal and the theory flows to an interacting Af = 2 
superconformal theory. Since the superconformal theory is nontrivial, one does not have a 
local effective Lagrangian description for those theories. The symmetry breaking pattern is 
obtained by the analysis at large /x >> A where there exists an effective description of the 
theory by integrating out the adjoint scalar <£>. The dynamical condensation of mesons breaks 
the USp(2N f ) flavor symmetry to U(N f ). 

2) The second group exists only when 2Nf > N — 4 5 with unbroken flavor symmetry 
USp(2Nf). Unlike the first group, the vacua can be in any of Higgs, Coulomb, or confining 
phases. They are given by the Special point in the SW-curve. In the Special point observed 
in [68], the gauge symmetry is enhanced to IR free SO(N) = SO(2Nf — N + 4) which is the 
Seiberg dual gauge group [75]. The full USp(2Nf) global symmetry remains unchanged since 
there are no meson condensates and no dynamical symmetry breaking occurs. These vacua are 
in the non-Abelian free magnetic phase. 

3.1.2 The strong coupling analysis 

When W(&) is very small compared with dynamical scale A, the M = 1 quantum theory 
can be considered as a small perturbation of a strongly coupled M = 2 gauge theory without 
W(&). In this region of parameters, we can use the Seiberg- Witten curve in which all the 
nontrivial dynamicses are already encoded. The M = 2 Seiberg- Witten curve derivation in the 
context of the matrix model was found in [38, 40]. Turning on the perturbed superpotential 
lifts most points on the Coulomb branch except some points (the Higgs roots) where a certain 
number of mutually local or nonlocal monopoles becomes massless. 

The strong coupling analysis has been done in [68, 64]. Let us recall that the curve of 

5 From this condition we see that when Nf = 0, it is satisfied only when N < 4. This explains why the Witten 
index of SO(N) without flavors is different for N < 4 because there are extra contributions from this second 
group. These special cases are also related to the fact that there is a factor x 4 (or x 2 ) in the Seiberg- Witten 
curve as we will see in the later examples. 
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SO(N) is given by 



[N/2] 



V 2 = II (t - ^ - AA^- N fh^ n (* - ml) 



3=1 



k=l 



where t 
reduced to 



x 



and e = for N odd and e = 1 for N even [76, 77]. For our massless case, it is 



y 2 = 



N/2 

JJ( f _^2 ) 2_ 4A 2(iV-2-iV / ) f JV /+ 2 j 
J'=l 

(JV-l)/2 

rj ( t _02 )2 _ 4A2 (^ 2 -^)^/+i ) 

i=i 



for N even, 



for N odd. 



By setting r 



6j's to vanish, the curve is 



V 



l2t 



\[N/2\-r 
3=1 



(3.2) 



which has a factor t with some power. The results in [68, 64] are that there are unlifted vacua 
under the mass deformation |/xTr$ 2 only when the power of factor t in the curve y 2 is some 
particular numbers, i.e., 



for N is even/odd or 



N c = (2N f -N + 4)/{2N f — + 3) 



(N f + 3)/(N f + l)/(N f + 2)/(N f + 2) 



for (N, Nf) = (e, e)/(e, o)/(o, e)/(o, o) where we denote e by an even number and o by an odd 
number. The first case [68, 64] is the character of the Special branch which in some sense 
corresponds to the baryonic branch of U (N) theory. The second case indicates the Chebyshev 
branch where the low energy effective M = 2 theory is a nontrivial conformal theory. 

The reason for the above conclusion is the following. For general r, it corresponds to the 
trivial super conformal theory classified as class 1 in [78]. To have the vacua, we must have 
enough massless monopoles which can happen if and only if r = N c = 2Nf — N + 4 when iV is 
even. It is easy to show that there are (N — Nf — 2) mutual local monopoles by noticing that 
from (3.2) (here we use the example where iV is even) 

y 2 = t 2N f -N+4 ( P *_ Nf _ 2 (t) - AA 2(N-2-N f ) t N-2-N f ^ 
_ ^2N f -N+4 ^j.N—2—Nf _ ^2(N-2-N f )y 



6 Notice that there is a factor (t — m|) instead of (t + mf,). The reason is as follows. The SU(N) case has a 
factor (x + TOfe) for det(a; + m) in the second term for diagonalized mass matrix m. For USp(2N) group, the 
masses are given by io-i ® diag(mi, m 2 , • • • , m,N f ) so det(x + m) contains a factor (x 2 + mf,). For SO(N) case, 
the masses are given by (2.3) so det(x + m) has a factor (x 2 — m\). 
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where to have a perfect square form in the bracket implying the maximal degeneracy of the 
Riemann surface, we have chosen P N _ Nf _ 2 (t) = t N ~ 2 ~ N f + A 2 ( N ~ 2 ~ N ^ [68]. This particular case 
is called the "Special branch" where the curve y 2 is a square form, which is a typical character 
of this branch. The Special branch, in fact, gives the second group of the solution analyzed in 
the weak coupling region. 

Besides the trivial superconformal fixed point, we can obtain a non-trivial superconfor- 
mal fixed point in IR by having a proper power of t in the curve, where mutually non- 
local monopoles are massless, and call it the "Chebyshev branch". To see these nontriv- 
ial fixed points, let us assume that both Nf and N are even. The curve is given by y 2 = 
t N f +2 (p 2 N _ Nf _ 2)/2 (t) - 4A2(A f -A f /-2)^ If we take the characteristic function as 



P(N-N f -2)/2(t) = 2( V AY N - N f-^T {N _ Nf _ 2) 
with r] 2< - N ~ N f~ 2 ^ = 1, we have the curve 7 

y 2 = t N f+ 2 A ^ A) 2(N-N f -2) 



'y/V 
2rjA / 



T 2 

(N—Nf—2) 



2r]A 



t N f+ 2 A ^ A y(N 



-N f -2) 



2r ] A / 



- 1 



(N—Nf—3) 



2r ] A / 



Notice that although rj 2< - N ' N f' 2 ^ = 1, the characteristic function P(N-N f -2)/2(t) is a function 
of t and 7] 2 , so we will have only the (N — Nf — 2) solution predicted from the weak coupling 
analysis. Notice also that except a factor t N i +2 , there are single roots at t — and t = A(r]A) 2 



from the factor 



\2r,\) 



U(N-N f -3) \ 2,) \ I ' 



(^jx)j so finally we have a factor t N f +3 for this case 8 . 
Notice that (Nf + 3) is an odd number which will be the property of the Chebyshev branch in 
the following discussions. Similar arguments can be made for other three cases. The Chebyshev 
branch discussed here will give the first group of the solutions in the weak coupling region. 



3.2 The factorized form of a hyperelliptic curve 

Now combining the proper factorization form of the SO(N) curve with the massless flavors 
described in the previous subsection with the proper factorization form of the U(Ni) curve 
with the flavors in [20], we can describe the general curve form. The basic idea is that first we 



need to have the proper prefactor (like t p for SO(N) part and (t 



m 



2\2r 



for the U(Ni) part at 



the r-th branch. More details about the power p will be presented in a later section.). After 

7 Hcre we have used a useful relation between the Chebyshev polynomials T£(x) — 1 = (x 2 — l)U 2 c _ 1 (x). 

They are defined as T~k{x) — cos(Kx) and Uk-i{x) — ■ 

8 The second kind of Chebyshev polynomial Utc{x) contains a factor x for odd degree K. 
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factorizing out these prefactors, we require the remaining curve to have the proper number of 
double roots and single roots, which will fix the form of factorization eventually. 

The key part of the above procedure is the number of single roots. The basic rule is the 
following. For every U (iVj) factor in the non-baryonic branch, we have two single roots while 
there is no single root in the baryonic branch. For the possible SO(N) factor, if it is at the 
Chebyshev branch, we have two single roots where one of them is at the origin; i.e., we have a 
factor t for that single root. However, if it is at the Special branch, there is no single root for 
the block SO(N). Adding the single roots together for all blocks, we obtain the final number 
of the single roots. If there are 2n single roots, we will have a factor F 2n (t) in the curve. 
Furthermore, as we will discuss soon, the polynomial F 2n {t) has some relationship with the 
deformed superpotential W'(x). One key point of all the above discussions is that the Special 
branch (or the Baryonic branch) will have one more double root than the Chebyshev branch 
(or the Non-baryonic branch) in the factorization form. This fact will be used repeatedly in 
the following study. 

To demonstrate our idea, we consider the following examples in which SO(2N c ) group is 
broken into the following three cases: 



where the generalization to multiple blocks will be trivial and straightforward 9 . Similar analysis 
corresponding to SO(2N c + 1) can be done with no difficulty. 

For the broken pattern SO(2N c ) — > SO(2N ) x U (Ni), by counting the number of the single 
roots, we derive the following four possible curves 



Curve (3.3) is for SO(2N ) at the Chebyshev branch and U(N 1 ) at the non-baryonic branch. 

The reason is that first we need have four single roots. Secondly, because SO(2N ) is at 

the Chebyshev branch, one of the four single roots must be at the origin t — and finally 

F 4 (t) = tF 3 (t). Curve (3.4) is for SO(2N ) at the Special branch and U(N 1 ) at the non-baryonic 

branch. Factor F 2 (t) with two single roots will record the information of U(l) C U(N 1 ). Curve 

(3.5) is for SO(2N ) at the Chebyshev branch and U(Ni) at the baryonic branch. Since 

SO(2N ) is at the Chebyshev branch we should have F 2 (t) = tFi(t). Finally curve (3.6) is 

9 For the breaking pattern SO(2N c ) — > U(N C ), the eigenvalues of $ are the same and are nonzero while for 
SO(2N c ) -> SO(2N c ), all of those are vanishing. 



SO(2N c ) -> SO(2N ) x U(N ± ), 



SO(2N c ) -> [/(Ay, 



SO(2N c ) -> SO(2N c ), 




tF 3 (t)H 2 Nr _ 2 (t), 
F^Hl^t), 
tF 1 (t)H 2 Nc _ 1 (t), 
H 2 Nc (t). 



(3.3) 
(3.4) 
(3.5) 
(3.6) 
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for SO(2N ) at the Special branch and U(N 1 ) at the baryonic branch, where no single root is 
required. Notice that the function H p (t) will have a proper number of (t — m 2 ) or t to count 
the prefactor for various branches. 

For the broken pattern SO(2N c ) — > U(N C ), we have 

y 2 = F 2 (t)H 2 Nc _ 1 (t), (3.7) 
V 2 = H 2 Nc (t), (3.8) 

where curve (3.7) is for U(N C ) at the non-baryonic branch and curve (3.8), for U(N C ) at the 
baryonic branch. Again, function H p (t) has a proper number of (t — m 2 ) to count the prefactor 
for various r-th branches. It is also interesting to see that the curve form (3.7) is identical to 
the (3.4), which as we will show later, provides the possibility for smooth transition between 
these two breaking patterns. 

For the broken pattern SO(2N c ) -> SO(2N c ), we have 

y 2 = tF 1 (t)H 2 Nc _ 1 (t), (3.9) 
y 2 = H 2 Nc (t), (3.10) 

where curve (3.9) is for SO(2N c ) at the Chebyshev branch and curve (3.10), for SO(2N c ) at 
the Special branch. The function H p (t) will have a proper number of factor t to count the 
prefactor required by Special or Chebyshev branch. It is noteworthy to notice that although 
both curves (3.5) and (3.9) appear to be identical, they can be distinguished by a factor H p (t) 
where different powers of t and (t — m 2 ) can arise. For example, there may be a factor (t — m 2 ) 
in H p (t) of the curve (3.5), but it does not exist in H p (t) of the curve (3.9). 

Before ending this subsection, let us give the rules for the counting of vacua. The total 
number of vacua is the product of the number of vacua of various blocks. For each block, there 
are four cases: U(Ni) without flavors, U(Ni) with flavors, SO(N) without flavors, and SO(N) 
with flavors. For U(Nj) without flavors, the counting is given by the Witten index w{N^) = iVj. 
For U (iVj) with Mj flavors, the result has been given in [20] and is summarized as follows: 



The number of vacua 



2Ni — Mi r< Mi/2, M : < N h 

Ni-Mi/2 r = M l /2, M % < N h 

2Ni-Mi r > M, - .V,, \I ; > X, ■ 1 

Ni-r r<Mi-Ni, M; > X, ■ 1. ' 

1 r = Ni — l(nonbaryonic), 

1 r = Mi- Ni, Ni (baryonic). 

For SO(N) without flavors, the counting is given by the Witten index w(N) = (N — 2) for 
iV > 5 and 4, 2, 1, 1, 1 for iV = 4, 3, 2, 1, (notice that in previous subsection we have given an 
explanation of the Witten index for N < 4). For SO(N) with M flavors, there are (N-M-2) 
vacua from the Chebyshev branch and one vacuum from the Special branch (if it exists). 
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4 Quartic superpotential with massive flavors 

In this section we will deal with and analyze the explicit examples with quartic tree level 
superpotential of degree 4 (that is, k — 1 from [2.2]) satisfying W'(±m) = 0, 

W'(x) =x(x 2 -m 2 ) (4.1) 

where we do not have a constant term in x because the superpotential W(x) is a function of 
x 2 with a constant term; therefore, after differentiating this with respect to x, superpotential 
W(x) will give rise to the above expression (4.1). In these examples the gauge group SO(N c ) 
can break into three cases: two blocks with 10 

SO(N c ) -> SO(N ) x u{n,), N c = N + 2N, 

under the semiclassical limit A — > which we will call the non- degenerated case and one block 
with 

SO(N c ) -> SO(N c ) or SO(N c ) -> U{[N~ c /2\) 

which we will call the degenerated case. Since for the sake of simplicity, we consider the case 
with equal masses of flavors, we can write the Seiberg-Witten curves for these gauge theories 
as follows: 

y 2 = P 2 Nc (x)-4x 4 A 4N *- 2N f- 4 (x 2 -m 2 ) N f for SO(2N c ), 
= P 2 Nc ( x )-4x 2 A m ^ 2N f- 2 (x 2 -m 2 ) N f for SO(2N c + l), 

where m in these curves is the same as the one in (4.1) 11 . On the r-th branch these Seiberg- 
Witten curves are factorized as follows, based on the discussion of previous section 3 and 
[73, 20], 

(x 2 -m 2 ) 2r [p| (iVc _ r) (x)-4x 4 A 4iVc - 2 ^- 4 (x 2 -m 2 )^- 2r ] for SO{2N c ), (4.2) 
( x 2 _ m 2yr [p* {Nc _ r) ( x ) - 4x 2 A 4N <- 2N f- 2 (x 2 - m 2 ) N f- 2r ] for SO(2N c + 1), (4.3) 

where the (x 2 — m 2 ) 2r is the prefactor for that particular branch. For the expressions in the 
bracket we need to count the proper number of single roots and double roots as given in previous 

10 As in [20], we use the notation for hat in U{N\) to denote a gauge theory with flavors charged under the 
U(Ni) group. 

11 For gauge theories with massless flavors we will study in next section in which the tree level superpotential 
(5.25) is different from (4.1) and all the flavors must be charged under SO (No), not U(N\). For degenerated 
case, SO(N c ) — ► SO(N c ), the flavors are charged under the factor SO(N c ). 
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section 3. Furthermore, there are two special r's which give the baryonic branch of U (N\) factor 

12 

r = iVi, or r = N f -N 1 . (4.4) 

If the r does not satisfy this condition, it is called a nonbaryonic branch. 

For convenience we list the proper factorization form of the curve in various situations. Note 
that the roots for W'(x) are x = and x = ±m from (4.1). For the non-baryonic branch of 
U(Ni) factor with some D5-branes wrapping around the origin x — 0, the curve corresponding 
to (3.3) is 

Pi {Nc -r) ~ 4* 4 A 4JVc - 2JV '- V - m 2 ) N f- 2r = x 2 H 2 2Nc _ 2r _ 4 (x)F 6 (x) for SO(2N c ), (4.5) 
Pm c - r) ~ 4x 2 A^- 2N f- 2 (x 2 - m 2 ) Nf - 2r = x 2 H 2 2Nc _ 2r _ 4 (x)F 6 (x) for SO(2N c + 1). (4.6) 

For the baryonic branch of U(Ni) factor with some D5-branes wrapping around the origin x = 
(or for the degenerated case where all D5-branes are wrapping around the origin x — 0, the 
curve has the same form: (3.9)), the curve corresponding to (3.5) is 

P 2 (Nc _ r) - Ax 4 A^- 2N f~\x 2 - m 2 ff- 2r = x 2 H 2 Nc _ 2r _ 2 (x)F 2 (x) for SO(2N c ), (4.7) 
P 2 (Nc _ r) - 4x 2 A^~ 2N f- 2 (x 2 - m 2 ) N f- 2r = x 2 H 2 2Nc ^ 2r _ 2 (x)F 2 (x) for SO(2N c + 1). (4.8) 

For the non-baryonic branch of U(Ni) factor (or for degenerated case with all D5-branes wrap- 
ping around the x = ±m, the curve is the same: (3.7)) the curve corresponding to (3.4) is 

Pi {Nc -r) - 4* 4 A 4JV - 2A ^ V - m 2 ff- 2r = H 2 2Nc _ 2r _ 2 (x)F 4 (x) for SO(2N c ), (4.9) 
P?(N c -r) ~ Ax 2 \^- 2N i~ 2 {x 2 - m 2 ) N f- 2r = H 2 2Nc _ 2r _ 2 (x)F 4 (x) for SO(2N c + 1). (4.10) 

Finally for the baryonic branch of U(N 1 ) factor (or for degenerated case with all D5-branes 
wrapping around the x = ±m, the curve takes the same form: (3.8) or (3.10)) the curve 
corresponding to (3.6) is 

^ e -r)-4x 4 A^- w /^-m 2 )^-* = J H^ e _ ar (x) for SO(2N c ), (4.11) 
P| ( ^_ r) -4x 2 A 4 ^- 2 ^- 2 (x 2 -m 2 )^- 2 ^// 2 V 2 ,W for SO(2N c+ l). (4.12) 

In fact, these curves have been given in section 3.2 where we used the fact that t = x 2 and 
r = 13 . 

Next we want to clarify the number of flavors Nf that we will study below. For the exponent 
of A in (4.2) and (4.3) to be positive, we will concentrate on S'0(2A^ C ) gauge theories with the 

12 This holds also for the degenerated case: SO(N c ) — > U([N c /2}). In this case, Ni becomes N c /2. 
13 Moreover, the curves corresponding to (3.7), (3.8), (3.9) and (3.10) can be written similarly and we will 
describe them in the following discussions when we need to explain each case. 
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condition N f < 2N C - 2 and SO(2N c + 1) gauge theories with the condition N f < 2N C - 1, 
which are the conditions for the theories to be asymptotically free. Thus taking into account 
of r = min (iV c , ^f^J , we will study the examples with the condition 

r < -J-. (4.13) 

Finally let us emphasize the comments on the relation between function Fq(x) appearing in 
(4.5) and (4.6) and W'(x) from (4.1). According to relation (A. 6), we have to take care of 
the second term of (A. 6) on the left hand side when the number of flavors Nf is greater than 
(2A C - 4) for SO(2N c ) gauge theories and (2A C - 3) for SO(2N c + 1) gauge theories. This 
is a kind of new phenomenon that did not appear in the pure case [19]. For example, on the 
non-baryonic branch of U(Ni) factor for 5*0(4) with Nf = 1, 5*0(5) with Nf = 2, 50(6) with 
Nf = 3 and 50(7) with Nf = 4 cases, we must be attentive to these modifications due to the 
presence of flavors. 

Now we are ready to deal with the explicit examples for SO(N c ) gauge theory with massive 
flavors where N c = 4, 5, 6, 7. The number of flavors Nf is restricted to Nf < N c — 2 and the 
index r satisfies (4.13). Therefore, for given number of colors N c , the quantities Nf and r are 
fixed. 

4.1 SO (4) case 

For this gauge theory we will discuss the number of flavors Nf = 0, 1 cases. There are three 
breaking patterns 50(4) -> 50(2) x 17(1), 50(4) -> 50(4), and 50(4) -> [7(2). 

4.1.1 N f = 

1. Non-degenerated case 

In [19], the case with Nf = was discussed. Since the factorization problem was trivial, 
the characteristic function P±{x) could be represented as 

P A (x) = x 2 (x 2 -v 2 ). 

From this P^x), we could obtain a tree level superpotential and a deformed function f2(x), 

W'(x) = x{x 2 - v 2 ), f 2 (x) = -Ax 2 A 4 . 

There is only one vacuum for a given W'(x). To satisfy the condition W'(±m) = discussed 
above, we should have v 2 = m 2 . Under the semiclassical limit A — > 0, the gauge group 50(4) 
breaks into 50(2) x U(l). Notice that the Witten index w(2) of 50(2) is one and therefore 
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the total number of vacua becomes one from the weak coupling analysis. We realize that 
the number of vacua from the strong coupling analysis coincides with the one from the weak 
coupling analysis. 

In fact, by factoring out x 2 common to W'(xf and f 2 (x), the full curve can be written as 

y 2 = x 2 (w\xf + f 2 (x)) = [(x 2 - v 2 ) 2 - 4A 4 ](x 2 ) 2 = F 4 (x)H 2 (x), 

which indicates that the 5*0(2) is at the Special branch where the power of x 2 is two (even) in 
the above curve and that U(l) at the nonbaryonic r = branch because the condition (4.4) is 
not satisfied. See also (4.9). 
2. Degenerated case 

The above calculations for nondegenerated case assume that every root of W'(x) has at 
least one D5-brane wrapping around it. However, there are the situations where some roots do 
not have wrapping D5-branes around them. For our example, there are two cases; one is that 
there is no D5-brane around zero (x = 0) and the other is that there is no D5-brane around 
x = ±m. For both cases, the curve corresponding to (3.7) should be factorized with factor 
F 4 (x) as 

y 2 = P 2 (x) - 4AV = F 4 (x)H 2 (x), 

where F 4 (x) can have a factor x 2 further depending on whether the SO(N) factor is at the 
Chebyshev branch or not. See also (3.9). Writing P 4 (x) = (re 4 — six 2 + s 2 ), H 2 {x) = (x 2 + a) 
and F 4 (x) = (x A + bx 2 + c), we have the following solutions 

& = 2(a^2A 2 ), c = a 2 , s 1 = -2a ± 2A 2 , s 2 = a 2 . 

If A — * 0, but a ^ 0, we get a symmetry breaking 50(4) — > U(2) where U(2) is at the 
nonbaryonic r = branch. If the classical limit goes to A, a — > where F 4 (x) contains the x 2 
factor, we get 50(4) — > 50(4). It is obvious that there is no smooth interpolation between 
these vacua because there exists a U(l) at the IR for the former symmetry breaking while there 
exists no U(l) for the latter symmetry breaking. 

To count the vacua, we need to determine the value of a. If there are D5-branes wrapping 
around the x = ±m, we get b = —2m 2 . Putting it back to the expression for b, we found two 
values for a (each sign gives one solution for a); thus, there exist two vacua for 50(4) — > U(2). 
If all D5-branes wrap around the origin x — 0, there are two cases to be distinguished. The first 
case is that the 50(4) factor is at Chebyshev vacua, so we need a factor x 2 in the curve which 
determines the value of c = (so a = too). There are two vacua corresponding to the ± sign. 
It is noteworthy that the curve is y 2 = (x 2 ) 3 (x 2 + b) with the power of x 2 being three. This is 
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the character of SO (A) without flavors at the Chebyshev branch. More detailed discussion of 
the power of t = x 2 can be found in the next section. The second case is that the SO (A) factor 
is at the Special vacua, so F±{x) should be a complete square form. This determines the value 
of a as a = A 2 if b = 2(a - 2A 2 ) or a = -A 2 if b = 2 (a + 2A 2 ). Although it seems that we have 
two solutions, in fact, the curve is the same for both cases and presumably we should count 
them as one. 

However, it is known that for SO (A) = SU{2) L x SU(2) R , we should have 2x2 = 4 
vacua. Among these four vacua, two of them have the relations S L = S R and S L + S R ^ 0, 
so they correspond to the above two solutions at the Chebyshev vacua. The other two have 
the relation Sl + S R — 0, but we found only one solution for the curve at the Special vacua. 
Our understanding is that since Sl + S R = for these two vacua, they correspond to the same 
point in the Af = 2 curve. Later we will meet same problem again and again, where from the 
curve we find only one solution for 5*0(4) at the Special branch, but it should be counted as 
two vacua. 

4.1.2 Nf — 1 

In this case, from the condition r < ^- we have only the r = branch. For the breaking 
pattern SO (A) — > SO (2) x U(l), the r = branch can be a non-baryonic or baryonic branch. 
On the other hand, for the breaking pattern SO(A) — > U(2), the r = is non-baryonic since 
r^Nf-N-L and r ^ Aq. 

• Baryonic r = branch 
Since we are considering a baryonic branch, we have an extra massless monopole (that is, one 
more double root) and from the relation (4.7) the following factorization problem arises, 

y 2 = p 2 ( x ) - Ax*A 2 (x 2 - m 2 ) = x 2 Hl(x)F 2 (x). 

To be consistent on both sides, the characteristic function P^x) on the left hand side must have 
an x 2 factor, i.e., P^{x) = x 2 (x 2 — a 2 ) due to the presence of an x 2 factor in the right hand side. 
In addition, we assume H 2 (x) = (x 2 — b 2 ). Assuming that 6^0 (the b = case will correspond 
to SO(A) — > SO(A) and will be discussed later), we must have F 2 (x) = x 2 , a 2 = [m 2 — A 2 ), and 
b 2 = (m 2 + A 2 ); thus, the characteristic function and F 2 (x) behave as 

p 4 (x) = x 2 [x 2 + A 2 - m 2 ) , F 2 {x) = x 2 . 

Under the semiclassical limit A — > 0, the gauge group SO (A) breaks into SO (2) x U(l). In this 
case, we have only one vacuum. This number matches the one obtained from the weak coupling 
analysis, because both SO (2) and U(l) with Nf — 1 theories have only one vacuum. Notice 
also that for this special case, the curve y 2 has a total square form with a factor x A implying 
that SO (2) is at Special branch. 
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• Non-baryonic r = branch 
1. Non-degenerated case 



In this case the factorization problem (4.5) becomes as 

P 2 (x) - 4x 4 A 2 (x 2 - m 2 ) = x 2 F 6 (x). 

This equation is easily solved if we assume that P±(x) = x 2 P 2 (x) = x 2 [x 2 — A) by recognizing 
that there is a factor x 2 on the right hand side and finally one obtains 

F 6 {x) = x 2 [(x 2 - A) 2 - 4AV - m 2 ) 

Taking into account of (A. 6) we have m 2 = A. Under the semiclassical limit A — > 0, the 
characteristic function goes to P±(x) = x 2 (x 2 — m 2 ), which means that the gauge group 5*0(4) 
breaks into 5*0(2) x U(l). We have only one vacuum, which matches the counting obtained by 
general analysis given in previous section (weak coupling analysis). The curve y 2 has a total 
square form with a factor x 4 implying that 50(2) is at Special branch. 
2. Degenerated case 

The curve should be factorized as, together with a term F 4 (x) on the right hand side as 
before, 

y 2 = Pl{ x ) - 4x 4 A 2 (x 2 - m 2 ) = F A {x)Hl{x). 

Parameterized by P^x) = (a; 4 — Six 2 + s 2 ), H 2 (x) = (x 2 + a) and F 4 (x) = (a; 4 + bx 2 + c) we 
have the following relations 

-2a - b - 4A 2 

51 = 2 ' 

4c + 4a (b - 4A 2 ) - (6 + 4A 2 ) 2 - 16A 2 m 2 

52 " 8 ' 

/ 32a 2 A 2 + (b + 4A 2 ) 3 + 16A 2 (b + 4A 2 ) m 2 + a (-2b 2 + 166A 2 + 32A 2 (3A 2 + m 2 ))\ 
° ~ ~{ 8a - 4 (6 + 4A 2 ) J' 

= a(4a 2 -4a(&-5A 2 ) + (5 + 4A 2 ) 2 ) 

+ ((-2a + bf + 8 (4a + b) A 2 + 16A 4 ) m 2 + 16A 2 m 4 . 

To see the limit where there is D5-brane wrapping around x = ±m, we put 14 b = —2m 2 — 4A 4 

14 To discuss the smooth transition, we take a to be free parameter while the mass m is determined by a, i.e., 
m will change as the a does. This is the philosophy used in [18, 20] (see, for example, the equation (6.30) in 
[20]) which is very convenient for this purpose. We can also use the direct method which will be demonstrated 
later by one simple example. To count the vacua, we really need to fix m and then solve a. 
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by the relationship of F A (x) and (x 2 — m 2 ) 15 and get 

= a 2 (a + 9A 2 ) + 3a (a + 4A 2 ) m 2 + (3a + 4A 2 ) m 4 + m 6 . 

There are two limits we can take: (1) A — > 0, but a — > constant. In this limit, all three solutions 
provide m 2 -> -a and it is 50(4) -> [7(2); (2) A -> 0, a -> which gives 50(4) -> 50(4). 
However, from the above equation by keeping the last term m 6 , we see that all of the solutions 
give m 2 — > in this limit. Since we require m 2 to be nonzero (that is, we are considering 
massive flavors), this limit is inconsistent with our assumption and should not be taken. In 
fact, as we have seen, the non-baryonic r = branch of 50(2) x U(l) theory will have a factor 
x A while for U{2) theory, the curve does not have any factor x 2 , so we should not expect any 
smooth transition between them. To count the number of vacua, we fix m and then solve a. 
There exist three solutions. They match the counting, (2iV - N f ) = (2 x 2 - 1) = 3, of [7(2) 
with one flavor. 

If all D5-branes wrap around the origin x — 0, there are two cases to be considered. The 
first one is when the extra x 2 factor occurs in F 4 (x): F^x) = x 2 (x 2 + d) or c = 0. In fact, it 
is the same curve of the baryonic r = branch, but with H 2 (x) = x 2 . There are two solutions 
with the expression P±{x) = x 2 (x 2 ± 2 A 2 ) by ± sign providing a breaking 50(4) — > 50(4) 
under the semiclassical limit A — > 0. These two vacua are the Chebyshev vacua of 50(4) gauge 
theory (notice that the counting becomes N — Nf — 2 = 4 — — 2 = 2) with an overall factor 
x 6 in the curve. The second one is that function F 4 (x) should be a complete square form. 
There are two solutions. One of them is Pa(x) = x 2 (x 2 — [m 2 — A 2 ]) which gives a breaking 
50(4) — > 50(2) x U(l) under the semiclassical limit A — > (it is the baryonic r = branch 
found above) and other, Pa(x) = + A 2 x 2 — m 2 A 2 which provides a breaking 50(4) — > 50(4) 
under the semiclassical limit A — > and 50(4) is at the Special branch. In other words, among 
the four vacua of 50(4), two of them are found in the Chebyshev vacua and two of them 
correspond to the same point in the Special vacua. 

Now we summarize what we have obtained in Table 1 by specifying the flavors Nf, symmetry 
breaking patterns, various branches, the exponent of t = x 2 in the curve, U(l) at the IR, the 
number of vacua, and the possibility of smooth transition. It turned out that the number of 
vacua is exactly the same as from the weak coupling analysis . 

4.2 SO(5) case 

For this gauge theory we will discuss the number of flavors Nf = 0, 1,2 cases. There are 

three breaking patterns 50(5) -> 50(3) x [7(1), 50(5) -> 50(5) and 50(5) -> [7(2). 

15 If we use t = x 2 , the curve is effectively the one for U(2) with three flavors where two of them are 
massless and one is massive. Then using the relation Fi(x) + 4A 2 x 2 ~ (x 2 — m 2 ) 2 + 0(x°), we can read off 
F 4 (x) =x a - (2m 2 + 4A 2 )x 2 + to 4 + d. 
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Table 1: The summary of the phase structure of the 50(4) gauge group. Here we use C/S 
for Chebyshev or Special branch for SO(Nj) factor and r NB /r B for the r-th non-baryonic or 
r-th baryonic branch. In this table, we list the power of t and the U(l) which is present in 
the nonbaryonic branch, at the IR. They are the indices to see whether two phases could have 
smooth transition. In the case of massive flavors, the power of t is 3 for SO(Ni) for iVj even 
and 1 for SO(Ni) for iVj odd at the Chebyshev branch. For massive flavors, only SO(Ni) where 
Ni = 2, 3, 4 can be at the Special branch and the power is 2, 0, for Ni = 2, 3, 4. For this table, 
we see that for both Nf = and Nf — 1, there are no two phases having the same factor t p 
where p is some number and the number ofU(l), so there exists no smooth transition between 
them. 

4.2.1 Nf — 

1. Non-degenerated case 

The theory with Nf = was discussed in [19]. Since the factorization problem was trivial, 
the characteristic function was given by, 

P 4 (x) =x 2 (x 2 -l 2 ). 

To satisfy the condition for W'(±m) = 0, we must have the relation, I 2 = m 2 . Under the 
semiclassical limit A — > 0, the gauge group 50(5) breaks into 50(3) x U(l) where 50(3) is at 
the Chebyshev branch with a factor x 2 in the curve while U(l) is at the nonbaryonic branch. 
However, since the Witten index w(3) = 2, we expect to find two solutions instead of one. In 
fact, another solution will be found in the degenerated case 16 . 

2. Degenerated case 



16 For the two vacua of 50(3), one is at the Chebyshev branch having factor t — x 2 and another, at the 
Special branch with factor t°. However, since the Special branch has extra double roots just like the one in the 
degenerated case, they have the same factorization form of the curve and for simplicity we use the degenerated 
case to include all of them. 
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For D5-branes wrapping around only one root of W'(x) (either the origin or x = ±m), the 
curve is 




V 2 = Pi it) ~ 4A 6 t = F 2 (t)H 2 (t), t = x 2 . 

First let us discuss a case in which there are D5-branes wrapping around ±m. Parameterizing 
^2(0 = {t 2 — Sit + s 2 ), F 2 (t) = (t — a) 2 + b and Hi(t) = (t + c), we found only one combination 
(si, s 2 , a, c) as a function of the parameter b (it is easy to see that b 7^ 0) 

b 2 4A 6 

b 2 4A 6 
a ~ 16A 6 ~ IT' °~ 16A 6 ' 

Now we can discuss the smooth interpolation. First using the relationship of F 2 (t) and (t — m 2 ) 
we get a = m 2 . If A -> 0, but 6 2 /A 6 = /3 ^ 0, a = m 2 -> (3/16 and P 2 (t) -> (i - ^) 2 
so the breaking pattern is 5*0(5) — > C/(2). If A — > 0, but A 6 /6 = 7 7^ 0, a = m 2 — > —47 
and P 2 (0 -> t(t + 47), the breaking pattern is 50(5) -> 50(3) x U{1) where 50(3) is at 
the Special branch. From this we see that U(2) is smoothly connected to 50(3) x U(l). It 
is noteworthy that at the Special point of 50(3) without flavor, the power of t is given by 
(2Nf — A^ + 3) = (0 — 3 + 3) = according to the argument in section 3.1.2. 17 Because of this, 
both 50(3) x U(l) and U(2) can have the same form of curve. 

Now let us count the number of vacua. To do so, we need to fix the mass m 2 and solve 
b. Using the relation a = m 2 , we find three solutions for b. It is obvious that two of them 
correspond to the limit b 2 /A 6 = (3^0 which gives the two vacua of U(2) and the third one, to 
the limit A 6 /6 = 7 7^ which gives the vacua of 50(3) x U(l). The third one is the missing 
one for 50(3) mentioned in the non-degenerated case. 

The above discussion shows the smooth transition clearly. However, to deepen our under- 
standing, we use another method by solving the b in terms of the m 2 . One of these three 
solutions will turn out 



-4(j) ^ 2 JT 



2 1 
2 (-) 3 (-9A 12 + 0V 18 (27A 6 - 4m 6 )) 3 . 



(-9A 12 + v^A 18 (27A 6 - 4m 6 
Let us consider the following expression in the above 

(-9A 12 + ^A 18 (27A 6 - 4m 6 )) 3 = A 3 (-9A 3 + V / 3 V / (27A 6 - 4m 6 ))" 

= 12^A 3 mc; 6 , w% = -1, A -> 

17 The proper power of t at the Chebyshev branch and Special branch will be summarized in the next section, 
where more details can be found. 
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where we have taken the limit A — * in the last line, but carefully kept the phase factor ujq. 
The phase factor comes from the expression 

(27A 6 - 4m 6 ) s = (4(m 2 - afj (m 2 - a 2 ) (m 2 - a 3 )) * 

where three roots aij are the order A 2 . If we start from very large m 2 , surround only one root 
once and go back to large m 2 , there will be a phase factor e - ^ = rj where if = 1. Thus 
depending on the path, in general we have 

b ~ -^=A 3 m[cu 6 r] k + c^VI- 

For general k, we find [u 6 r] k + VI ^ and b ~ A 3 . This case gives the breaking pattern 
5*0(5) — > C/(2) in the limit A — > 0. However, when A; = 1, we have {ooqI]) 2 = if = — 1 (where 
we have used cj| = 77) and [cj 6 r] fc + ujQ 1 r}~ k ] = 0. This case will give rise to b ~ A 6 and the 
breaking pattern 50(5) — > 50(3) x £7(1). From this analysis we see that starting from k = 
and surrounding one root once, we make a smooth transition from one phase to another phase. 

The above analysis demonstrates explicitly the origin of the smooth transition among the 
phases in the parameter space of deformed superpotential. However, in most cases it is hard 
to make such a detailed analysis. It is more convenient to use b as a free parameter while m 2 is 
determined by b when our focus is just the smooth transition. Later we mainly use this method. 
To count the number of vacua, we really need to fix m 2 and find the number of solutions of 
b. Although we cannot solve this exactly most of the time, the counting can be done easily by 
observing the degree of the equation and various limits we take. 

Finally we consider a case in which all D5-branes wrap around the origin x — 0. In this 
case, we need to have a factor t in F 2 (t), or a 2 + b = using the above results 18 . There 
are three solutions b = — 4A 4 , 2(1 ± iy / 3)A 4 which give three vacua with the breaking pattern 
50(5) -> 50(5) where 50(5) is at the Chebyshev branch. The counting goes like (N - 2) = 3 
which is consistent with the one from both the strong and weak coupling analyses. 

4.2.2 N f = 1 

In this case, there is only one nondegenerated breaking pattern 50(5) — > 50(3) x U(l) 
where we can have both non-baryonic and baryonic r = branch due to the fact r = Nf — Ni 
and r 7^ Ni. Besides that, we have two degenerated breaking patterns 50(5) — > 50(5) 
and 50(5) — > U(2) where only the non-baryonic r = branch exists because the relations 
r ^ N f - Ni and r ^ N x hold. 

18 Notice that for this case, a, b are fixed by this equation a 2 + b = 0, so there is no free parameter in the curve 
to be adjusted. This explains why the above smooth transition, where b is a free parameter, does not include 
this phase. 
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• Baryonic r = branch 

The curve should be factorized together with a function F 2 (x) as 



y 2 = p 2 ( x ) - Ax 2 A\x 2 - m 2 ) = x 2 H 2 (x)F 2 (x). 
Setting H 2 (x) = (x 2 — b), F 2 (x) = x 2 — a and Pi(x) = x 2 (x 2 — si), we get 



, 2A 4 m 2 4A 4 m 2 , 4 2 a4 4 , 

si = 6 — — , a = — — , - 6 4 + b 3 m 2 + A 4 m 4 = 0. 4.14 

tr tr 

From (4.14) we see two limits: (1) If A — > and 6 — > constant, by keeping the first and second 
terms we find m 2 — > 6 and the gauge group is broken to 5*0(5) — > 50(3) x £7(1) where 50(3) 
is at the Chebyshev branch and U(l), the baryonic r = branch; (2) If A — > 0, but 6 ~ A 4 / 3 , 
by keeping the second and the third terms we find m 2 — > 6 3 /A 4 ~ finite and the gauge group is 
broken to 50(5) — > 50(5). Thus, we see a smooth transition between these two phases. This 
is also consistent with the indices. Both phases have the same t l and zero U(l) factor at IR. 

Now we can count the vacua by solving b for fixed m 2 . Clearly there are four solutions. 
Keeping the first and second terms we get one solution b ~ —m 2 which has the classical limit 
50(3) x £7(1) where 50(3) is at the Chebyshev branch. Keeping the second and third terms 
we get three solutions b ~ A 4 / 3 which give the three vacua of 50(5) and are consistent with 
the counting — 2 = 3 from the weak coupling analysis. 

To show clearly the rightness of the above result, we use the m 2 to solve b directly. One of 
these four solutions will be 

A = -9A 4 m 8 - v / 3A 4 mV256A 4 + 27m 4 , 
As m 4 42h^A 4 m 4 

2j — 




Let us consider the classical limit A — > 0. First we have A < A 4 . Using this factor and 
considering the phase factor as we did for the 50(5) with Nf = 0, we have v^S ~ rf^ with 
i] 2 = 1. Using this, we finally have 

(l-77-^2(l-77)), L0 2 = 1 



4 

771 



4 



with two phase factors 77,0; depending on the path of m 2 . For 77 = 1 we have b ~ at the 
classical limit. For 77 = —1, u — 1 again b ~ 0. However, for 77 = —1, uo = — 1 we have b ~ m 2 . 
The limit of b ~ gives the breaking 50(5) — > 50(5) while the limit b ~ m 2 gives the breaking 
50(5) — > 50(3) x U(l). These calculations confirm our conclusions derived by the simpler 
method. 
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For another vacuum of SO (3) x U(l) where SO (3) is at the Special branch, we just need 
to choose Pa(x) = x 2 (x 2 — m 2 ) + A 4 so that y 2 = (x 2 (x 2 — m 2 ) + A 4 ) 2 is a total square form. 
• Non-baryonic r = branch 

1. Non-degenerated case 
The curve should be 

y 2 = p 2 ( x ) - Ax 2 A 4 (x 2 - m 2 ) = x 2 F 6 (x). 

It is easy to solve and it turns out P^(x) = x 2 (x 2 — s) and F 6 (x) = x 2 [(x 2 — s) 2 — 4A A (x 2 — m 2 )}. 
To determine the parameter s, by noticing that for Nf = 1, F 6 (x) = W'(x) 2 + ax 2 + b, we get 
s = m 2 . Here we found only one solution for 5*0(5) — > SO (3) x U(l) where SO (3) is at the 
Chebyshev branch with a factor x 2 in the curve. Another one will be given from the Special 
vacua of SO (3) in the degenerated case. 

2. Degenerated case 

For the degenerated case where D5-branes wrap around only one root, the curve should be the 
form with a factor F±(x) 

y 2 = p 2 ( x ) - Ax 2 A 4 (x 2 -m 2 ) = F 4 (x)Hl(x). 

Parameterizing as P±{x) = (x 4 — s±x 2 + s 2 ), H 2 (x) = (x 2 + a), and F\{x) = (x 2 — 6) 2 + c, we 
found the following solutions 

-2ab + c + 4A 4 4A 4 (a-b + m 4 ) 

si = -a + b, s 2 = , c 



2 a + b 

= (a 2 + am 2 )((a + 6) 2 -4A 4 )-A 4 m 4 . 

Using the relationship between F±{x) and (x 2 — m 2 ), we get b = m 2 . Putting it into the last 
equation above we get 

= (a + m 2 fa - 4A 4 a(a + m 2 ) - A 4 m 4 . (4.15) 

This equation (4.15) has three solutions of m 2 . If A — > 0, but a — > constant, keeping the 
first and third terms we have all three m 2 going to —a and (a + m) 2 ~ A 4 / 3 . This limit gives 
SO (5) — > U(2). If A — > and a — > but a ~ A 4 , keeping the first and third terms we get one 
m 2 -> A 4 /a ^ (another two solutions of m 2 go to 0). This limit gives SO (5) -> 50(3) x 17(1) 
where 50(3) is at the Special branch. From this discussion we see that at least one solution of 
m 2 gives a smooth transition between these two phases. 

To count the vacua, notice that for fixed m 2 we have four solutions a from (4.15). Keeping 
the first and third terms we see that there are three solutions a — > — m 2 and one solution a — > 0; 
i.e., there are three vacua for 50(5) — > U (2) by counting 2A^ — Nf = A — 1 = 3 and one vacuum 
for 50(5) -> 50(3) x [7(1) where 50(3) is at the Special branch. 
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4.2.3 N f = 2 

In this case we have the following situations. For 5*0(5) — > 5*0(3) x U(l), we can have the 
r = non-baryonic branch with two vacua and the r = 1 baryonic branch with two vacua. 
For 50(5) — > 50(5), there are three vacua, while for 50(5) — > U(2) we can have the r = 
non-baryonic branch with two vacua, the r = baryonic branch with one vacuum, and the 
r = 1 non-baryonic branch with one vacuum. We will study these vacuum structures in detail. 

• Baryonic r = 1 branch 
As previously discussed, on r = 1 branch, the Seiberg-Witten curve is factorized as follows: 



The solution for this equation is easily obtained. Pi(x) = x 2 because this characteristic function 
should have an x 2 factor in order for both sides to be consistent with each other. Thus, we have 
P 4 (x) = x 2 (x 2 — m 2 ), which means that the gauge group 50(5) breaks into 50(3) x U(l). For 
fixed to, there is only one vacuum which comes from the Chebyshev point with a factor x 2 in 
the curve. To get another vacuum, we must choose P2(x) = (x 2 + A 2 ) , so that 



This vacuum comes from the Special point. Thus, we get all of two vacua for r = 1 baryonic 
branch. 

• Non-baryonic r = 1 branch 

In this case, we need use the relationship of F^{x) = P2(x) 2 — 4x 2 A 2 and (x 2 — m 2 ) 2 so that 
p2(x) = (x 2 — to 2 ) . There is only one solution which gives the one vacuum of 50(5) — > U(2) 
at the r = 1 non-baryonic branch. 

• Non-baryonic r = branch 
1. Non-degenerated case 

In this case, from (4.6) we have the following factorization problem, 



To satisfy this equation we must have x 2 factor in P±{x). Thus, assuming that P±{x) = x 2 [x 2 — 
a 2 ), we can obtain F 6 (x) as follows: 




In addition to this factor, we have the following factorization problems, 



Pl(x) -Ax 2 K 2 = x 2 F 2 {x). 



P 2 (x)-4x 2 A 2 = (x 2 -A 2 ) 2 . 



P!(x)-4x 2 A 2 (x 2 -m 2 ) 2 



x 2 F 6 (x). 



F 6 (x) = x 2 {x 2 - a 2 ) 2 - 4A 2 (x 2 - to 2 ) 2 . 
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At this point we want to determine W'(x) from F 6 (x). As already discussed above, we should 
take into account the second term of (A. 6) 19 . 

F 6 (x) + 4AV (x 2 - m 2 ) = W\xf + 0(x 2 ). 

From this equation, we can read off W'(x) = x 2 (x 2 — a 2 ). Thus, according to the condition 
W'(±m) = 0, we identify a 2 with m 2 . Finally the characteristic function becomes 

P A (x) = x 2 (x 2 - m 2 ) , 

which means that the gauge group 5*0(5) breaks into 50(3) x £7(1) where 50(3) is at the 
Chebyshev branch with a factor x 2 in the curve. For fixed m, there is only one vacuum. We 
will find another vacuum coming from the degenerated case. 

2. Degenerated case 
Now let us consider the degenerated case with the curve (note the factor F 4 (x)), 

y 2 = p^ x ) - 4x 2 A 2 (x 2 - m 2 ) 2 = F 4 (x)H 2 (x). 

Writing Pa(x) = (x A — s\x 2 + S2), H 2 (x) = (x 2 + a) and F 4 (x) = x A + bx 2 + c we found 
-2a - b - 4A 2 



S2 = 



4c + 4a {b - 4A 2 ) - (6 + 4A 2 )' - 32A 2 m 2 



I 3 2a2 A 2 - 2a(& -l 2 A 2 y + 4A 2 ) + (, + 4A 2 ) 3 + ^ ^ + fc + ^ ^ + 

a - I - 2A 2 

'4a 2 - 4a (b - 5A 2 ) + (b + 4A 2 )^ + 24aA 2 m 2 + 4A 2 m 4 . 
Putting b = —2m 2 — 4A 2 we can solve 20 

= 4a 2 (a + 9A 2 )+8a(a + 3A 2 )m 2 + 4(a + A 2 )m 4 ^m 2 = 

At the limit A — > and a — > constant, we get both m 2 — > —a and 5*0(5) — > £7(2). At the limit 

A -> 0, but -> iA + aA 3 , we have m 2 ^ and 50(5) -> 50(3) x 17(1) where 50(3) is 

at the Special branch. Thus, there is smooth transition between these two phases. To count 

the vacua, notice that for fixed m 2 , we have three solutions for a. Setting A = 0, the equation 

becomes 4a(a + m 2 ) 2 = 0. Thus we have two vacua for U(2) and one vacuum for 50(3) x U(l). 

19 In this case, we are considering the tree level superpotential (4.1), thus gm+2 is 1. In addition, since the 
equation (A. 6) is the result for SO(2N c ) theory, we have to extend the result to SO(2N c + 1) theories by 
replacing ^4^-2^-4 in thc sccon d tcrm to p^N c -2N f -2^ 

20 In this case, c = — 4aA 2 + 4A 2 m 2 + m 4 is very simple compared with the above messy expression. 
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Next we discuss the case in which all D5-branes wrap around the origin x = 0; thus, -F 4 (x) 
should have a factor x 2 . Writing Pi(t) = t(t — si), Hi(t) = t + a and F 2 (t) = t(t + b) (where 
we have used t = x 2 for simplicity), we get s 1 = -a - 2 A 2 + 2A ^™ 4 , b = — 4A ^ 4 and 

= a 2 (a + A 2 ) - 2aA 2 m 2 + A 2 m 4 , or m 2 = ±«-^~ + «• 

There is only one limit a ~ A 2 / 3 to get finite m 2 . For fixed m 2 , there are three solutions of a 
which give the three wanted vacua of 5*0(5) — > SO (5) where SO (5) is at the Chebyshev branch 
with a factor x 2 = t in the curve. 

Finally we want to get the r = baryonic branch of U(2). To achieve this, we require y 2 to 
be square. See also (4.12). There are two solutions. One of them is Pi{i) — (t — m 2 )(t + A 2 ) 
(for simplicity we have used t = x 2 ) which gives the r = 1 baryonic branch of 5*0(3) x U(l). 
Another is p2(t) — (t — m 2 ) 2 + A 2 t which is the one we want. It is easy to see that the 
latter solution gives y 2 = [(£ — m 2 ) 2 — A 2 t] 2 ^ if 1 (t) 2 F 1 (t) 2 , which explains why we could not 
get it in the previous paragraph. In other words, for the baryonic branch, the relationship 
F 4 (x) + 4AV = (^) 2 + O(x ) fails. 

Now we summarize the results in Table 2 by specifying the flavors Nf, symmetry breaking 
patterns, various branches, the exponent of t = x 2 in the curve, U(l) at the IR, the number 
of vacua, and the possibility of smooth transition. It turned out that the number of vacua is 
exactly the same as from the weak coupling analysis. 

4.3 SO(6) case 

The next example is 50(6) gauge theory, which is more interesting. For this gauge theory 
we will discuss the number of flavors Nf = 0, 1,2,3 cases. There are three breaking patterns 
50(6) -> 50(4) x 17(1), 50(6) -> 50(2) x [7(2), 50(6) -> 50(6), and 50(6) -> U{3). 

4.3.1 Nf — 

1. Non-degenerated case 

In [19] the theory with Nf = was discussed intensively. For the breaking pattern 50(6) — > 
50(2) x U(2) where 50(2) is at the Special branch, there are two vacua which are confining 
vacua. These confining vacua are constructed from the Coulomb branch with the breaking 
pattern 50(4) — > 50(2) x U(l) by the multiplication map. In addition to these vacua, there 
are two vacua with the breaking pattern 50(6) — > 50(4) x U(l). These two vacua come from 
the Chebyshev branch of 50(4) factor, and we will find another two vacua from the Special 
branch of 50(4) in the degenerated case. 

2. Degenerated case 
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Group 


Branch 


Power of t(= x 2 ) 


U(l) 


Number of vacua 


Connection 





SO(3) x 17(1) 


(C,0 NB ) 


t l 


1 


1 




(S,0 NB ) 


t° 


1 


1 


A 


SO(5) 


(C) 


t L 





3 




U(2) 


(Onb) 


t u 


1 


2 


A 


1 


SO(3) x U(l) 


(C,0 NB ) 


t 1 


1 


1 




(S,0 NB ) 




1 


1 


B 


(C,0 B ) 


¥ 





1 


D 


(S,0 B ) 







1 




SO(5) 


(C) 


¥ 





3 


D 


U(2) 


(Onb) 


t° 


1 


3 


B 


2 


SO{3) x U(l) 


(C, NB ) 


t 1 


1 


1 




(S,0 NB ) 




1 


1 


C 


(C, 1 B ) 


¥ 





1 




(S, 1 B ) 







1 




SO(5) 


(C) 


¥ 





3 




U(2) 


(Onb) 


t° 


1 


2 


C 


(Ob) 


t° 





1 




(Inb) 




1 


1 





Table 2: The summary of the phase structure of SO(5) gauge group. There are three pairs of 
smoothly connected phases indicated by the letters A,B,C. That is, the two branches denoted 
by A, B, or C are smoothly connected with each other. For Nf = 2, since we can have r = 0, 1, 
the r is also an index for phases and should be preserved under the smooth transition. One 
application is that although for Nf = 2, SO(3) x U(l) at (S, 1 B ) and U(2) at (0 B ) have the 
same t° and zero U(l), they can not be smoothly connected because the index r is different. 
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To see all the vacua including the one with the breaking pattern 5*0(6) — > U(3), it is 
necessary to consider the degenerated case. The factorization problem becomes 

P^(x) - Ax 4 A 8 = Hj(x)F 4 (x). 

After solving this factorization problem, we obtain the following solutions, 

j-j / \ 6 3G 4 + A 8 4 / 2 3G 6 \ 2 G(G 4 -A 8 ) 2 



/ r' 4 + a 8 \ 2 

= K~^H + (4-16) 



where e 2 = —1. Note that this use of e is slightly different from previous one. If we remember 
the relation F±(x) = W ^ + d, we can see the relation, 

A = e ^^T- i wh ere W'(x) = x (x 2 - A) . (4.17) 

In this case we can take two semiclassical limits with A — > 0: 

1. G — > with fixed = t>: Under these limits, since the characteristic function goes 
to - ► (^ 2 ~ V T i the gauge group 50(6) breaks into U(3). From the condition (4.17), 

3 1. 

A = e^j, we obtain three values for eG = (— AA 4 ) 3 . Thus substituting (4.16), we can see three 
vacua, which agree with the number obtained by weak coupling analysis, because pure U(3) 
theory has three vacua. 

2. G — > with fixed = w: Under this limit since the characteristic function behaves 
as Pq(x) — > x 4 (x 2 — w), the gauge group 50(6) breaks into 50(4) x U(l). From condition 
(4.17), A = e^, we obtain two values for G. These two vacua of 50(4) x U(l) are the 
missing two vacua mentioned in the non-degenerated case. Adding all the contributions from 
the degenerated case and nondegenerated case, the result agrees with the one from the weak 
coupling analysis, because the Witten index for pure 50(4) gauge theory has four vacua. 

Notice that from the above analysis we get a smooth transition between U(3) and 50(4) x 
£7(1) at the Special branch. The reason is that the power of t = x 2 of 50(4) at the Special 
branch is (27V^ — + 4) = 0, which is the same as the case of U(3). Similar phenomena have 
been observed in the previous subsection of the 50(5) gauge group. 

Finally, to get the breaking 50(6) — > 50(6) (see also (3.9)), we require a x 2 factor in F^(x) 
which can be obtained when G 4 — A 8 = 0. There are four solutions for G. Combining the e, 
it seems we will get eight solutions. However, careful study shows that there are, in fact, only 
four solutions which match the counting (N — 2) = 4 of the Witten index for 50(6) gauge 
theory where 50(6) is at the Chebyshev branch. 
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In this example, there is only the r = branch. For the breaking pattern 5*0(6) — > 
50(4) x U(l), this r = branch can be non-baryonic or baryonic due to the fact r = Ni 
while for the breaking pattern 50(6) — > 50(2) x U(2), this branch is non-baryonic because 
both relations r ^ Nf — Ni and r ^ Ni hold in this case. In addition, for the breaking 
pattern 50(6) — > U(3) it is a non-baryonic branch because r 7^ Nf — N\ and r 7^ iVi. For 
50(6) — > 50(6), only the Chebyshev branch exists with four vacua. 

• Non-baryonic r = branch 

1. Non-degenerated case 



Since this branch is non-baryonic we have the following factorization problem, 

P 2 (x) - 4x 4 A 6 (x 2 - m 2 ) = x 2 H 2 (x)F 6 (x). 

To satisfy this equation, the characteristic function must have an x 2 factor, Pq{x) = x 2 P^(x). 
Then one of H 2 (x) and F Q (x) on the right hand side must have a factor x 2 . If the factor x 2 is 
contained in the function Fq(x), it belongs, in fact, to the degenerated case, so we consider the 
case H 2 (x) = x 2 . Parameterizing Pq{x) = x 2 (x A — Six 2 + s 2 ) and F 6 (x) = x 6 + ax A + bx 2 + c, 
we found the solutions: 

ct a 2 
s 1 = — , s 2 = =F2iA 3 m, b = — =f AiA 3 m, c = -4A 6 =f 2iaA 3 m. 

Obviously the classical limit is 50(6) — > 50(4) x U(l). By the relationship of W'(x) and F 6 (x) 
we find a = —2m 2 , so there are two vacua of 50(4) x U(l) which come from the Chebyshev 
branch of the 50(4) factor where there is a factor x 6 in the curve. 
2. Degenerated case 



Next we go to the degenerated case. The factorization problem of this case is given by 

P 2 (x) - Ax 4 A e (x 2 -m 2 )= Hj(x)F 4 (x). 

Parameterizing P 3 (t) = t 3 - sit 2 + s 2 t + s 3 , H 2 (t) = t 2 - at + b, F 2 {t) = t 2 - 2m 2 t + m 4 + c 
where we have used the relationship between F 2 (t) and (t — m 2 ) 2 to simplify the calculation 21 . 
The solution is 

c 4aA 6 2 4 

— h am — m , 

4 c 



21 The logic of the calculation is following. First we factorize the curve with some parameters. Then we 
establish the relationship between these parameters and to 2 . This is the method used in 5*0(4) and SO(5). 
However, we can use the relationship of to 2 and ^4(2;) to solve some parameters firstly, then to solve the 
factorization form. The latter method is more convenient for our purpose. 
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a 



= 



2c (4cA 6 (c 3 - 64A 12 ) + c 2 (c 3 + 64A 12 ) m 2 + 8A 6 (3c 3 + 64A 12 ) m 4 ) 
(c 3 + 64A 12 ) (c 3 - 64A 12 + 32cA 6 m 2 ) 

((c 3 - 64A 12 ) 2 - 1024c 2 A 12 m 4 - 1024cA 12 m 8 ) . 



(4.18) 



From the last equation in (4.18), we see that at A — > we must have c — > 0. Therefore the 
solution of m 2 is roughly 



The first limit is that c ~ A 2 - 4 , so that all four solutions provide m 2 ^ and a ~ 2m 2 and 
b ~ m 4 . This limit gives 5*0(6) — > [7(3). The second limit is that c ~ aA 12 , so that all four 
solutions give m 2 ^ 0, a ~ ^A 6 ~ 0, 6 ~ 0. This limit gives 50(6) -> 50(4) x [7(T) where 
50(4) is at the Special branch. Thus, we have a smooth transition between these two phases. 

To count the vacua, for fixed m 2 there are six solutions of c. Keeping the first and third 
terms of (4.18) we have five c ~ A 2 - 4 , which is consistent with the counting (2N—Nf) = 6—1 =5 
of U(3). Keeping the second and third terms of (4.18) we have one solution of c ~ A 12 , which 
gives the one vacuum of 50(4) x U(l) where 50(4) is at the Special branch 22 . 

Now we consider other special cases. The first one is that b = so that H 2 (t) = t(t — a). 
There are three solutions: a = —uj 3 A 2 + m 2 , c = — Auo^ 1 A 4 with cjf = —1. They are the three 
vacua of 50(6) -> 50(2) x [7(2) by counting number (2N - N f ) = 4 - 1 = 3. 

The next one is that F 2 (t) = t(t + a) and H 2 (t) = t(t + b). The solutions are given by 



From the (4.19) we see that a — > in the limit A — > 0. There are two limits we can take. If 
a ~ 4A 3 / 2 , we have m 2 ~ —1, b ~ which gives 50(6) — > 50(6). If a ~ A 3 , we have m 2 ~ —16, 
6 ~ 16 which gives 50(6) -> 50(4) x Z7(l) where 50(4) is at Chebyshev branch and [7(1) is 
at the r = baryonic branch. 

From these considerations we see that there is a smooth transition between 50(6) — > 50(6) 
and 50(6) — > 50(4) x [/(l). Indices for these two phases are the same: t 3 and non IR U(l). 
It is interesting to compare them with the baryonic r = branch of 50(5) with Nf — 1, where 
we do not find a smooth transition between 50(5) — > 50(5) and 50(5) — > 50(3) x C/(l). For 
50(5), these two phases come from different solutions of m 2 while for 50(6) there is only one 
solution for m 2 . This different behavior of m 2 as a function of parameters explains the different 
phase structure, but the physics behinds these calculations is still unclear. 

To count the vacua at fixed m 2 , keeping the first and the third terms of the second equation 

(4.19) we find four solutions a ~ m 1 / 2 A 3 / 2 which give the limit 50(6) — > 50(6). Keeping 

22 Notice that although the Special branch of SO (A) should have two vacua, but the counting from the curve 
gives always one. We have met this situation many times already. 



(c 3 - 64A 12 ) 2 
1024cA 12 




a 16A 6 



a 6 + 4096A 12 + a 2 A 6 (128a + 256m 2 ) = 0. 



(4.19) 
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the second and the third terms we find two solutions a ~ ±4iA 3 /m and b ~ — m 2 . These 
two solutions give the breaking pattern 5*0(6) — > 5*0(4) x U(l) where 50(4) is at Chebyshev 
branch and U(l) is at the r = baryonic branch. 

• Baryonic r = branch 
Since this branch is baryonic, from the relation (4.7) we have the following factorization prob- 
lem, 



In fact, for the 50(4) at the Chebyshev branch, we need to have H 2 (t) = t(t + b). This is 
exactly the same problem discussed above. 

To get the missing vacua, we just need to choose P${t) = t 3 — m 2 t 2 + A 6 , so that y 2 = 
(t 3 -m 2 t 2 - A 6 ) 2 . This gives the r = baryonic branch of 50(6) -> 50(4) x [7(1) where 50(4) 
factor is at the Special branch. 

4.3.3 N f = 2 

In this case, we have two branches r = 0, 1. For the breaking pattern 50(6) — > 50(4) x U (1), 
the r = 1 branch is baryonic because r = Nf — Ni = N± while the r = branch is non-baryonic. 
For the breaking pattern 50(6) — > 50(2) x U(2) the r = 1 branch is non-baryonic while the 
r = branch can be non-baryonic or baryonic because r = Nf — Ni . Moreover, for the breaking 
pattern 50(6) — > U(3) both branches of r = 0, 1 are non-baryonic. For the breaking pattern, 
50(6) — > 50(6), there is only Chebyshev branch. 

• Baryonic r = 1 branch 
The Seiberg-Witten curve is factorized as, 



In addition to this factor, we have the following factorization problem with a factor F 2 (x), 



To satisfy this constraint, the characteristic function P^{x) must have x 2 factor, i.e. P±{x) = 
x 2 (x 2 — a 2 ). In addition, we assume that H 2 (x) = (x 2 — b 2 ). After solving the factorization 
problem we obtain the following results, 



where r\ is 2-nd roots of unity. Thus taking into account the (x 2 — m 2 ) factor, we obtain 



P 2 (t) - 4t 2 A 6 (t - m 2 ) = tHl^F^t). 




(4.20) 



y 2 = Pl( x ) - 4x 4 A 4 = x 2 H 2 2 {x)F 2 {x). 



P 4 (x) = x 2 (x 2 - 2r]A 2 ), F 2 (x) = {x 2 - 4r^A 2 ), 



P 6 (x) = x 2 (x 2 - m 2 )(x 2 - 2r]A 2 ). 
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Under the semiclassical limit, A — > 0, the characteristic function becomes Pq(x) — > x 4 (x 2 — m 2 ), 
which means that the gauge group 5*0(6) breaks into 5*0(4) x U(l). Since there is rj we have 
two vacua for fixed m which come from the Chebyshev branch of 50(4) factor. 
In addition to this solution, we have other solutions from the following case, 

Pt(x)-4x*A 4 = Ht(x). 

If we parameterize H 4 (x) = (x 2 — a)(x 2 — b), we have two solutions a — > —A 2 ,b — > A 2 or 
a — > A 2 , b — > —A 2 . Both cases lead to the conclusion that P 4 (x) becomes P 4 (x) = x 4 + A 4 . 
As above, these solutions also exhibit the breaking pattern 50(6) — > 50(4) x U(l). However, 
these come from the Special branch of 50(4). Therefore, we have four solutions from the 
strong coupling analysis. This number agrees with the one from the previous analysis of the 
weak coupling. Although 50(4) gauge theory has four vacua, U(l) with Nf = 2 gauge theory 
has only one vacuum. 

• Non-baryonic r = 1 branch 

1. Non-degenerated case 

As in the previous case, we can factorize the Seiberg- Witten curve as (4.20). In this case, 
since the branch is non-baryonic, we have different factorization problem, 

y 2 = Pl{x)-Ax A \ 4 = x 2 F,{x) 

in which the factor (x 2 — m 2 ) 2 has been factorized out. To solve this equation, we assume 
that P 4 (x) = x 2 (x 2 — a 2 ). After inserting this relation, we find that Fq(x) must have an x 2 
factor. However, this implies the degenerated case. Therefore, we have no new solution for this 
non-degenerated case. 

2. Degenerated case 

Next we will consider the degenerated case. The factorization problem is given by 

Pl{x) - 4a; 4 A 4 = H 2 (x)F 4 (x). 

Solving this factorization problem, we obtain two kinds of solutions. One is given by 

P 4 (x) = (x 2 - aj 2 + 2r]A 2 x 2 , 

F 4 (x) = [x 2 -a + 2r]A 2 ) 2 + 4ar]A 2 -4A\ 

where r\ is the 2-nd root of unity. Taking into account that W'(±m) = 0, we have one constraint, 

m 2 = a + 2r]A 2 . 
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In this case, we can take only one semiclassical limit, A — > 0. Under the limit, the characteristic 
function goes to P%(x) — > (x 2 — m 2 ) 3 , which means that the gauge group 5*0(6) breaks into 
U(3). Thus, we have two vacua from the strong coupling approach. This number of vacua 
agrees with the weak coupling analysis given in (3.11). The U(3) theory with Nf — 2 and r = 1 
has two vacua. 

The other solutions are derived from H 2 (x) = x 2 . Inserting P 4 (x) = x 2 (x 2 — a 2 ) into the 
factorization problem, we obtain 

F 4 (x) = (x 2 - a 2 ) 2 - 4A 4 = (^^j 2 + d. 

The condition W'(±m) = leads to the solution a 2 = m 2 . After all, we obtain the trivial 
solution, 

P 6 (x) = x 2 [x 2 -m 2 ) 2 . 

The gauge group SO (6) breaks into 5*0(2) x U(2). Thus, we obtained only one vacuum from 
the strong coupling analysis which comes from the Special branch of 50(2). This number 
matches the one from the weak coupling analysis. The U(2) with Nf = 2 gauge theory has 
only one vacuum from (3.11). Since 50(2) theory has only one vacuum, the total number of 
vacua is only one, which agrees with the results of the strong coupling analysis. 

• Non-baryonic r = branch 

1. Non-degenerated case 
Next we consider a non-baryonic r = branch. From (4.5) we have the following factorization 
problem, 

P 2 (x) - 4x 4 A\x 2 - m 2 ) 2 = x 2 H 2 (x)F 6 (x). (4.21) 

Assuming that P§(x) = x 2 P 4 (x), we can rewrite (4.21) as follows: 

x 2 (P 4 (x) - 2 A 2 (a; 2 - m 2 )) (P 4 {x) + 2 A 2 (a; 2 - m 2 )) = Hl{x)F & {x). 

If we include the additional factor x 2 in a function F 6 (x) it goes to the degenerated case. Here 
we want to discuss the non-degenerated case so we consider the H 2 (x) = x 2 case. The solution 
is given by 



Pdx) = x 2 



x 4 + Ax 2 - 2r]A 2 m 2 

2 , /*4 , a2 2\ ( 2 , A-2 V A 2 



F 6 (x) = x 2 (x 2 + A) -4(A 4 + r ? A 2 m 2 ) [x 2 + 



m 2 



m 2 + r]A 2 

where rj is the 2-nd root of unity. Taking into account that W'(±m) = we have the condition 
m 2 = —A. Under the semiclassical limit, A — > 0, the characteristic function goes to Pq(x) — > 
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x 4 (x 2 — m 2 ), which means that the gauge group 5*0(6) breaks into SO (A) x U(l). These vacua 
come from the Chebyshev branch of the 5*0 (4) factor. 
2. Degenerated case 

Now we consider the degenerated r = branch. Using t = x 2 , the curve should be factorized 

as 

P|(t) - 4A 4 t 2 (t - m 2 ) 2 = H 2 (t)F 2 (t). 

We have two cases to be discussed for this r = branch. The first case is H 2 (t) = tHi(t) which 
constrains to Ps(t) = tP 2 (t). Then we have the following solution 

P 2 (t) - 2 V A 2 (t - m 2 ) = H 2 (t) = (t- a) 2 , 

P 2 (t) + 2r]A 2 {t - m 2 ) = (t - a) 2 + 4r]A 2 {t - m 2 ), 

with rf — 1. There are two situations we need to consider for this kind of solution. The first 
is that there are D5-branes wrapping around the x = ±m; then we have 

a = m — 2rjA 2 , 

by the relationship of F 2 (t) and (t — m 2 ) 2 . It gives 50(6) — >■ 50(2) x U(2) at the classical limit 
where 50(2) is at the Special branch. There are two vacua corresponding to rj — ±1, which is 
consistent with the counting that U (2) with two flavors at r = non-baryonic branch has two 
vacua. The second situation is that F 2 (t) has a further factor t which leads to 

a 2 - Ar]A 2 m 2 = . 

It gives 50(6) — > 50(6) at the classical limits. There are four solutions, which is consistent 
with the counting that 50(6) at the Chebyshev branch has a Witten index (N — 2) = 4. Notice 
that because of the second situation, F 2 (t) has a factor t, and the first situation is not smoothly 
connected to the second situation. 

Having considered the case that H 2 (t) = tHi(t), we move to the case that there is no factor 
t in the curve y 2 . Then Ps(t) — 2r]A 2 t(t — m 2 ) has only one double root for every rj, or 

P 3 (t) + 2A 2 t(t - m 2 ) = (t-a-b) 2 (t + c), 

P 3 (t) - 2A 2 t(t - m 2 ) = (t- a + b) 2 (t + d), 

where 

c 

d 




= 26 + 2A 



, a{b + 2A 2 ) A 2 m 2 



= -26-2A 2 



+ 

o o 
a(b + 2A 2 ) A 2 m 2 
b +_ 6~' 
6 2 (6 + A 2 ) + aA 2 (m 2 -a). 
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Using the relationship of F 2 (t) = (t + c)(t + d) and (t 

2ab{b + 2A 2 ) - 26A 2 m 2 
&2 = 



m 2 ), we get 



2m 2 , 



from which we can solve 



a = 



(b + A 2 )m 2 
6 + 2A 2 



c = 2(6 + A 2 ) - m 2 , 



d = -2(6 + A 2 ) - m 2 



and 



(& + A 2 ) (b 2 + 



(6 + 2A 2 ) 2 



) 



= 0. 



(4.22) 



Forgetting about the factor (& + A 2 ) we found that there is only one limit b ~ A to get finite 
m 2 . In this limit m 4 ~ —1, a ~ m 2 , c,d~ —m 2 , so it gives SO (6) — > £7(3). For fixed m 2 , there 
are four solutions of b which give the four vacua of U (3) by counting {2N — Nf) = 6 — 2 = 4. 

However, where is the breaking SO (6) — > 5*0(4) x U(l)7 It comes from (b + A 2 ) = 0, so 
that a = and c — d — — m 2 . In fact, the curve has the same form as the curve of the r = 1 
branch. Similar phenomena have been observed in [20]; for example, the U(2) with Nf = 2. 
Since this non-baryonic r = branch is coincident with the baryonic r = 1 branch in this 
particular example, we do not have a smooth transition from U(3) to 5*0(4) x U(l). It is very 
interesting to recall that there is smooth transition for Nf = 0, 1. 

• Baryonic r = branch 

For the baryonic branch, the curve should be 



Parameterizing H^(t) = t 3 — at 2 + bt + c, there are two solutions. The first one is a = 2m 2 , 
b = m 4 — A 4 , and c = which gives r = baryonic branch of 50(2) x U(2). The 50(2) is at 
the Special branch with a factor t 2 in the curve. The second solution is a = m 2 , b = —A 4 , and 
c = — A 4 m 2 which gives the r = 1 baryonic branch of 50(4) x U(l) discussed before. 

4.3.4 N f = 3 

In this case, some of the branches can be described by the addition map we have discussed in 
previous section. 

• Non-baryonic r = 1 branch 

In this case, the curve is simplified as 



where the term in the bracket is, in fact, the curve of 50(4) with Nf — 1. This is just the 
familiar application of addition maps. Using the corresponding result of 50(4), we find three 



P 2 (t) - 4A 4 t 2 (t - m 2 ) 2 = Hl(t). 



y 2 = ( x 2 - m 2 ) 2 p 2 ( x ) - Ax 4 A 2 (x 2 - m 2 ) , 
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vacua for 50(6) -> U(3) and one vacua for SO(Q) -> 50(2) x £7(2). There is no smooth 
transition in this case. 

• Baryonic r = 1 branch 

Again, we can use the addition map to reduce the problem to the r = baryonic branch 
of SO (A) with Nf — 1. From there, we get one vacuum for 50(6) -> 50(2) x £7(2), two 
vacua for 50(6) — > 50(4) x £7(1) where 50(4) is at the Chebyshev branch, and two vacua for 
50(6) -> 50(4) x [7(1) where 50(4) is at the Special branch. 

• Non-baryonic r = branch 
1. Non-degenerated case 

In this case, we have the following factorization problem, 

Pl(x) - 4x 4 A 2 (x 2 - m 2 ) 3 = x 2 Hl{x)F & {x). 

To satisfy this equation, we must have a x 2 factor in Pq(x) = x^P^x). Let us assume that 

p A (x) = x 4 + Ax 2 + B, H 2 {x) = x 2 - C, 
F 6 (x) = [x(x 2 - D)f + G (x 2 + F) . 

After solving the factorization problem, we have three kinds of solutions. However, two of them 
have an x 2 factor in F 6 (x), which means that this becomes the degenerated case. Therefore, in 
the non-degenerated case, we have only the solution below. 
The third solution is given by 

P 6 (x) = x 2 [x A + Ax 2 - 2^Am 3 ] , = 4 (AA 2 - A 4 + 3 Am 2 - ir?Am 3 ) , 
m 3 (A - SirjAm) D = _ A + 2 A 2 



iArjA — ir/A 3 + 3ii]Am 2 + m? ' 

Thus taking into account the relation W'(x) and F 6 (x), we have m 2 = —A. Under the semi- 
classical limit, the characteristic function behaves as Pq(x) = x 4 (x 2 — m 2 ), which means that 
the classical group 50(6) breaks into 50(4) x U(i). Thus, we have obtained two vacua from 
the analysis of strong coupling, which comes from the Chebyshev branch of the 50(4) factor. 
2. Degenerated case 

For this case, the curve should be factorized as 

P 2 (t) - AA 2 t 2 (t - m 2 f = F 2 (t)H 2 (t). 

There are two cases we should discuss. The first case is that H 2 (t) = tHi(t). After dividing out 
the t 2 , at two sides we parameterize P 2 (t) = t 2 — Sit + s 2 , Hi(t) = (t — a), and F 2 (t) = t 2 + bt + c. 
Solving Si,s 2 ,c in terms of a and b we get the following constraint 

A 2 [a - m 2 ) 3 (4a 2 + 4a (b - 5A 2 ) + (b + 4A 2 ) 2 + 36A 2 m 2 ) = 0. 
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At the classical limit, we have b = —2a 2 . If there are D5-branes wrapping around the x = ±m, 
we will have b ~ —2m 2 , thus giving SO (6) — > 5*0(2) x U (2). For given 6, there are two solutions 
for a which are consistent with the counting of U{2) with Nj = 3 at the r = non-baryonic 
branch. If there are no D5-branes wrapping around the x = ±m, we need to have an extra t 
factor in F 2 (t)] i.e., c = 0. There are four solutions having a, b — > which give the four vacua 
of SO (6) -> 50(6). 

Now we move to the second case in which there is no factor t in H 2 (t). This case will 
give the r = non-baryonic branch of 50(6) — > C/(3) and the r = non-baryonic branch 
of 50(6) -> 50(4) x £7(1) where 50(4) is at the Special branch. Although the indices for 
these two phases are the same, it is not clear whether they are smoothly connected or not. To 
demonstrate the smooth transition, we must use the analytic expression, which is too complex 
to be solved. To count the number of vacua we can use a numerical method. Since we are 
unable to solve it, we will not discuss if any further. 

• Baryonic r = branch 

This can happen only for 50(6) — > U(3) and the corresponding curve should be factorized 

as 

P 2 (t) - AA 2 t 2 (t - m 2 ) 3 = H 2 (t). 

There are three solutions for H 3 (t) = t 3 — at 2 + bt — c. The first one is a = 2m 2 + A 2 , 
b = m 2 (m 2 + A 2 ), and c = which gives the r = 1 baryonic branch of 50(2) x U{2). The 
second one is a = m 2 + A 2 , b = 2A 2 m 2 , and c = A 2 m 4 which gives the r = 1 baryonic branch of 
50(4) x U(l). The third one is a = 3m 2 + A 2 , b = 3m 4 , and c = m 6 which gives the baryonic 
r = branch of U (3) we are looking for. 

The results are summarized in Table 3 by specifying the flavors Nf, symmetry breaking 
patterns, various branches, the exponent of t = x 2 in the curve, U(l) at the IR, the number 
of vacua, and the possibility of a smooth transition. It turns out that the number of vacua is 
exactly the same as from the weak coupling analysis. 

For the 50(7) case, we refer to version one in the hep-th archive for detail. In this paper 
we only summarize the results in Table 4 and Table 5 by specifying the flavors Nf, symmetry 
breaking patterns, various branches, the exponent of t = x 2 in the curve, U(l) at the IR, the 
number of vacua, and the possibility of a smooth transition. It turned out that the number of 
vacua is exactly the same as from the weak coupling analysis. 

5 Quartic superpotential with massless flavors 

Thus far we have discussed the phase structures of massive flavors. In this section we 
will focus on massless flavors. Classically, under the breaking pattern SO(N c ) — > SO(N ) x 
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N, 


Group 


Branch 


Power of £(= x 2 ) 


[7(1) 


Number of vacua 


Connection 





SO(2) x 17(2) 


(5,0^) 


t 2 


l 


2 




50(4) x 17(1) 


(C,0 NB ) 


t 3 


1 


2 




(5, Otvb) 


t u 


1 


2 


A 


50(6) 


(C) 


t 3 





4 




[7(3) 


(Oatb) 


m 

t u 


1 


3 


A 


1 


S0(2) x 17(2) 


(5, Oats) 


t 2 


1 


3 




50(4) x 17(1) 


(C.Ojvb) 


* 3 


1 


2 




(5,0^) 


t u 


1 


2 


B 


(C,0 B ) 


t 3 





2 


C 


(5,O b ) 


rn 

t u 





2 




50(6) 


(0) 


t 3 





4 


c 


[7(3) 


(Onb) 


t° 


1 


5 


B 


2 


50(2) x [7(2) 


(5, Ijvb) 


t 2 


1 


1 




(5,0 B ) 


t 2 





1 




(S,0 NB ) 


* 2 


1 


2 




50(4) x [7(1) 


(C, 1b) 


t 3 





2 




(5, 1b) 


t u 





2 




(C, NB ) 


t 3 


1 


2 




(5, Oivs) 




1 


2 




50(6) 


(0) 


"73 

t 3 





4 




[7(3) 


(Inb) 


*° 


1 


2 




(0 NB ) 


t u 


1 


4 




3 


50(2) x [7(2) 


(5, Ijvb) 


t 2 


1 


1 




(5, 1b) 


t 2 





1 




(S,0 NB ) 


t* 


1 


2 




50(4) x [7(1) 


(C, 1b) 


t 3 





2 




(5, 1b) 


t° 





2 




(O, Ojvb) 


t 3 


1 


2 




(5, Otvb) 


t u 


1 


2 


* 


50(6) 


(C) 


* 3 





4 




[7(3) 


(Ijvb) 


t° 


1 


3 




(Ojvb) 


t u 


1 


3 


* 


(Ob) 


t U 





1 





Table 3: The summary of phase structure of 50(6) gauge group. Here again we use capital 
letters to indicate the three phases having smooth transitions. The * means that it is not clear 
whether they are smoothly connected or not. 
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N, 


Group 


branch 


Power of t(— x 2 ) 


[7(1) 


Number of vacua 


Connection 





50(3) x £7(2) 


(S,0 NB ) 


t° 


l 


2 


A 






(C,0 NB ) 


t 1 


1 


2 


B 




50(5) x 17(1) 


(C,0 NB ) 


t 1 


1 


3 


B 




50(7) 


(O) 


, 1 

t l 





5 






(7(3) 


/ pi \ 


t 


i 
i 


Q 

o 


A 

A 


1 


50(3) x [7(2) 


(0, WB ) 


. i 


1 


3 


c 






(5, O^b) 


t u 


1 


3 


D 




50(5) x 17(1) 


(0,(W> 


t 1 


1 


3 


C 






(C,0 B ) 


t 1 





3 


N 




50(7) 


(C) 


t 1 





5 


N 




U(3) 


(Onb) 


t° 


1 


5 


D 


2 


50(3) x 17(2) 


(O, Iatb) 


t 1 


1 


1 








(5, Ijvb) 


t° 


1 


1 


E 






(O, O^vb) 


t l 


1 


2 


F 






(5,0^) 


t° 


1 


2 


G 






(5,0b) 


t v 





1 








(O,0b) 


t L 





1 







50(5) x 17(1) 


(0,1b) 


t 1 





3 








(O,0iVB) 


t 1 


1 


3 


F 




50(7) 


(C) 


t 1 





5 







[7(3) 


(Inb) 


t° 


1 


2 


E 






(Onb) 


t" 


1 


4 


G 


3 


50(3) x [7(2) 


(C,0 NB ) 


t 1 


1 


1 


H 






(S,0 NB ) 


t° 


1 


1 


I 






(O, Iatb) 


t L 


1 


1 








(5, Ijvb) 


t u 


1 


1 


J 






(O, 1b) 


t l 





1 








(5, 1b) 


t» 





1 






50(5) x [7(1) 


(O, Ojvb) 


t 1 


1 


3 


H 






(0,1b) 


t l 





3 






50(7) 


(C) 


t 1 





5 






[7(3) 


(Oatb) 


t° 


1 


3 


I 






(Ijvb) 


t° 


1 


3 


J 






(Ob) 


t° 





1 





Table 4: Summary of the phase structures of 50(7) gauge group. Here again we use capital 
letters in the last column to indicate the various phases having a smooth transition. 
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A T 

N f 


Group 


Branch 


rower or t(— x ) 


U(l) 


AT v, r 

IN umber of vacua 


Connection 


4 


£0(3) x U(2) 


(C 1 , Ijvb) 


t 1 


1 


1 








^, J- ATS J 




1 


1 


K 






(C, NB ) 




1 


2 


L 






(S, NB ) 




1 


2 


M 






\y, ^b) 


/l 


n 
u 


i 
i 








IA Z B) 


I 


n 
u 


i 
i 






SO(o) x (7(1) 


(C, 1 B ) 


J.1 





3 








(C,0 NB ) 


■i 

r 


-I 

1 


3 


T 

L 




SO(7) 


(C) 







5 






U(3) 


(Inb) 


i° 


1 


2 


K 






(Onb) 


t u 


1 


3 


M 






(2>nb) 


t u 


1 


1 








(1b) 


t u 





1 





Table 5: Continued. Summary of the phase structures of SO (7) gauge group. 

Y\^ =1 U(Nj), these massless flavors will be charged under SO (No) since W'(x)\ x= ± m =o = is 
always true in our considerations 23 . 

To demonstrate the new features coming from the massless flavors and compare them with 
the weak and strong coupling analyses, let us start with the following quadratic superpotential 
first 

= -/iTr$ 2 . 

The curve is given by 

where e = for N c odd and 1 for iV c even and Nf < N c — 2 for asymptotic free theory. For the 
non-baryonic branch, the curve should be factorized as 

y 2 = F A (x)H\ Nc/2] _ 2 (x) (5.23) 

with Fi(x) = x 2 F2(x) where the appearance of x 2 in this case is characteristic of the non- 
baryonic branch. See for example, (4.7) and (4.8). Furthermore, we have the following rela- 
tionship 

F 2 (x) = W'(x) 2 + d = x 2 + d . 
23 For degenerated cases, we have SO(N c ) — > SO(N c ) and SO(N c ) — > U([N c /2]) for quartic superpotential. 
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Because of the special form of F 4 (x), we do not have the r-th branch in the U(N C ) case. 

To show it more clearly, let us consider the example of SO(2N c ) gauge theory with 2M 
flavors. If we use t = x 2 , the curve is given by 

y 2 = f[(t- $f - 4A 2(2JV c -2-2M) t 2 + 2M > 
3=1 

The non-baryonic vacua should be factorized as 

y 2 = F 2 (t)H 2 Nc _ 1 (t), F 2 (t)=t(t + d) . 
Due to the factor t of F 2 (t), we should have, for example, 4>n c = and the curve becomes 

= t(t + d)H 2 (t) 



y 2 = t 2 



N c -1 

]J(t- (f) 2 ) 2 - 4A 2 ( 2iVc - 2 - 2M )t 2+2 ( M - 1 ) 

3=1 



which means that the function HN c -i(t) must have a factor t. Cancelling out factor t 2 at both 
sides, we get 

JVc-l 

y\ = f[( t - 2 ) 2 - 4A 2(2iV c -2-2M) t 2 + 2(M-l) = f(f + d ) H ^ t) 
3=1 

which is the non-baryonic curve of SO(2(N c — 1)) with 2(M — 1) flavors. Since we require that 
the number of flavors satisfies 2M < 2N C — 2 for asymptotic free theory, the above operation 
can be iterated at most (M + 1) steps and the factorization problem is reduced to 

y 2 = t 2 ^[P 2 Nc _ {M+1) (t) - AA 2 ^- 2 - 2 ^] = t(t + a)^ M ^H%_ {M+1) _M 

The presence of factor t 2 ( M+1 ) indicates that no matter what factor t 2r with r < M + 1 we start 
with, we always end up with r = M + 1 branch. This particular branch is nothing but the 
Chebyshev point emphasized in the weak and strong coupling analyses given in section 3.1.2. 
It is easy to check that for other N c and Nf, the factorized curves are also at the Chebyshev 
point. For the Special branch, we have already written down the form of curves in the strong 
coupling analysis where the basic feature of this branch is that the curve is a complete square. 

Having obtained the quadratic superpotential, we can progress to the general superpotential 
with degree 2(n + 1) where k = n. Again, for example, the SO(2N c ) with 2M flavors at the 
general non-baryonic branch should have the factorization, by generalizing (5.23) 

V 2 = ft(* - ^?) 2 - 4A 2 ^- 2 ~ 2M h 2 + 2M = tF 2n+1 {t)H 2 NcHn+1) {t). 

3=1 

Similar to the above derivation, factor t tells us that at least one of 0j's is vanishing, so that 
there is a factor t 2 on the left hand side. Then either ifjv c -(n+i)(i) or F 2n +i(t) must be divided 
by factor t. 
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If t is present in the F 2n+1 (t), we should have F 2n+ i(t = x 2 ) = x 2 ( W ^ ) 2 + f n ~i(x 2 ) and 
the curve becomes 

y 2 = F 2n (t)[tH N _ {n+1) (t)} 2 (5.24) 

with, (for example see (A. 9)), 

F 2n (t = X 2 ) = P^) 2 + fn^X 2 ) . 

The factorization form (5.24) is, in fact, the expression for the Special branch because there is 
no overall factor t in the curve. Geometrically, the function F 2n+1 {t) divided by t means that the 
cycle around the origin t = on the reduced Riemann surface is degenerated and the remaining 
non-trivial cycles are around U(N) factors. In other words, the non-baryonic solution found 
in this case comes from the Special root discussed in the weak and strong coupling analyses. 

Instead, if t is present in the function H N _^ n+1 ^(t), then we can cancel the t 2 factor on both 
sides of the curve and reduce the problem to SO(2(N c — 1)) with 2(M — 1) massless flavors. 
Iterating this procedure, we can reach the Chebyshev point as in the quadratic superpotential 
if it does not stop at the above Special point. 

In summary, the above analysis shows that the factorization of curves will stop only at 
either the Chebyshev point or the Special point, which is consistent with the results from the 
weak and strong coupling analyses. 

So considered, we can now focus on the following quartic superpotential W(Q>) = |Tr$ 4 — 
^Tr$ 2 , so that 

W'(x) = x(x 2 - a 2 ) . (5.25) 

We also consider only the case Nf < N c — 2 for asymptotic free theory. For convenience 
of the following discussions, we would like to summarize the essential points of the form of 
Seiberg-Witten curves: 

a) . At the Chebyshev point, the power of t = x 2 in the curve is (Nf + 3) for (Nf,N) = 
(even, even), (N f + 1) for (Nf, N) = (even, odd), (N f + 2) for (Nf, N) = (odd, even), and 
(Nf + 2) for (Nf,N) = (odd, odd). Notice that for all four cases, the power is always an odd 
number which is consistent with the factorization form tH 2 (t); i.e., except one factor of t, the 
other power of t can be considered as double roots. The number is given by (N — Nf — 2). 

b) . At the Special point which exists only when N c = 2Nf — N + 4 > 0, the characters of the 
curve are: (1) When N f < N - 2, the power of t is (2N f — N + A) for iV even and (2N f — N + 3) 
for N odd. When Nf > N — 2, the power of t is iV for N even or (N — 1) for iV odd. Notice 
that the power is always even, so it can be considered to belong into the double root in the 
factorization form of curve; (2) the number of the double root increases by one compared with 
that at the Chebyshev point. 
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c) . From (a) and (b) the simplest way to distinguish between the Chebyshev branch and 
the Special branch is to see whether the power of t is odd or even. 

d) . Because of point (b), the Special branch has a similar form of factorization as one of 
degenerated cases and will be discussed together in the following discussion. 

5.1 50(4) case 

Because of Nf < N c — 2 for asymptotic free theory, we only have two choices: Nf = and 
Nf — 1. The Nf = case has been discussed in section 4 under massive flavors where we only 
need to replace m by a. However, the Nf = 1 will give new results. There are three classical 
limits S0{4) -> S0{2) x 17(1), SO(4) -> S0{4), and SO {A) -> 17(2). For SO(2) with one 
flavor, there exists only the Special branch which can also be seen by taking into consideration 
the counting of the Chebyshev branch, N — Nf — 2 = 2— 1 — 2<0, which does not exist. For 
5*0(4) with one flavor, we can have both the Chebyshev branch and the Special branch because 
it is possible that 2N f - N + 4 = 2 > and N - N f - 2 = 1 > 0. 

1 Non-degenerated case 
The curve should be 

y 2 = Pt(x)-4x 6 A 2 =x 2 F (i (x) 

by taking the m — > limit into (3.3). Using P±(x) = x 2 [x 2 — A), we get Fq(x) = x 2 [(x 2 — A) 2 — 
Ax 2 A 2 ]. From the relationship between Fq(x) and W'(x) we find that A = a 2 . In the classical 
limit A -> 0, we get SO(A) -> SO(2) x U(l). Notice that there is a x 4 = t 2 factor in the y 2 
indicating that the solution found comes from the Special root instead of the Chebyshev point. 
2. Degenerated case and Special branch 

As previously noted, since the Special branch has one extra double root and the degenerated 
case has one root of W'(x) without wrapping D5-branes, the starting factorization forms of the 
curve are the same and we should consider them simultaneously. The curve should be 

y 2 = p 2 ( x ) - 4a; 6 A 2 = F A (x)H 2 (x) 

or 

y 2 = [t 2 - sit + s 2 ) 2 -4t 3 A 2 = (t 2 + at + b)(t + c) 2 , t = x 2 . (5.26) 

Given the a there are three solutions for (s u s 2 , 6, c). The first is ( ~ a ~ 4A2 , 0, ( ~°~ 4A2 ) 2 , 0), where 
c = gives a factor t 2 in the curve. Notice that there is no smooth transition found in this 
solution. 

To determine the classical limit, we need to determine the behavior of a. There are two 
cases to be analyzed. The first case is that there are some D5-branes wrapping around x = ±a. 
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Group 
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Table 6: The phase structure of SO (A) with one massless flavor. There are no smooth transitions 
between these vacua because there is no vacua which has the same number of the power of t 
and the number ofU(l). 

Then we get a = —(2a 2 + 4A 2 ), which produces the classical limit 5*0(4) — > 50(2) x U(l). 
In fact, this solution is identical to the one given in the previous paragraph where we claimed 
that it belongs to the Special branch. There is only one solution which is consistent with the 
counting. The second case is that all D5-branes wrap around the origin x = 0, so curve (5.26) 
describes the 50(4) — > 50(4) classically. With one flavor, we can have the Special branch with 
a factor t 2 by (2Nf — N + 4) = 2 or the Chebyshev branch with t 3 where Nf + 2 = 3. To get 
the Special branch we must set a = — 2A 2 so that (t 2 + at + b) — (t — A 2 ) 2 . There is only one 
solution. To get the Chebyshev branch, we need b = so that a = — 4A 2 . Notice the counting 
of vacua matches the formula (N — Nf — 2 = 4—1 — 2 = 1). 

The second and third cases have a somewhat complex expression as a function of a and 
neither b or c is zero for general a 

-2a + A 2 =f 3AV-2a + A 2 

51 = 2 ' 

a 2 - 19aA 2 + 16A 4 ± 5aAV-2a + A 2 =f 16A 3 V-2a + A 2 

52 = 1 ' 

b = — + 2A 4 =f 2AV-2a + A 2 + a(-3A 2 ± AV-2a + A 2 ), 

a - 5A 2 ± 3AV-2a + A 2 
C = 2 • 

Under the limit A -> 0, (s 1 ,s 2 ,b,c) -> (-a, X' T' f ) which S ives 50(4) -> U{2) with two 
vacua coming from the U(2). Again by the relationship between F^{x) and W'(x)/x we get 
a = —(2a 2 + 4A 2 ). Notice that we cannot take the limit A — > 0, a — > because the curve 
form (5.26) has no factor t before the limit is reached, whereas if there are D5-branes wrapping 
around origin x — 0, the starting curve must have a factor t. Because of this inconsistency, we 
can not take that limit. Similar observations have been presented in [20]. 
The phase structure has been summarized in Table 6. 
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5.2 50(5) case 

For the SO (5) gauge group, Nf can be zero, one, or two. There are three breaking patterns 
SO (5) -> 50(3) x 17(1), 50(5) -> 17(2), and 50(5) -> 50(5). The N f = case has been 
discussed in section 4 under massive flavors. For the Nf = 2 case, notice that the curve is 

y 2 = P!(x) - 4AV 

which is exactly the same form as the curve of 50(4) with iV/ = 1. Thus, the results are 
the same as 50(4) with one flavor. In fact, it is an example of the addition map previously 
discussed in Appendix B. For Nf — 1, the curve 

y 2 = P 4 2 (x)-4AV 

is the same as that of 50(4) without any flavors, so the results can be applied here. However, 
now we need to distinguish between Special vacua and Chebyshev vacua for the case of massless 
flavors. Because of this, we will repeat the 50(5) example, but not discuss 50(7) at all 
because everything in 50(7) can be reduced to the results of 50(6). Additionally, we will use 
a somewhat different method to redo the calculations. It will be interesting to compare these 
two methods. 

1. Non-degenerated case 
The curve should be 

y 2 = p 2 ( x ) - 4:r 4 A 4 = x 2 F 6 (x). 

Using Pa{x) = x 2 (x 2 — u) and F e (x) = W'(x) 2 + ax 2 + b, we found that b = 0, u = a 2 , and 
a = —(a 4 — 4A 4 ), which gives one vacuum of 50(5) — > 50(3) x U(l). Again, the fact that 
y 2 = x^F^x) indicates that it is the Special vacuum. The reason that the Chebyshev branch 
of 50(3) x U(l) does not exist can be seen from the counting A^ — Nf — 2 = 3— 1 — 2 = for 
the Chebyshev branch. 

2. Degenerated case and Special branch 
The curve should be 

y 2 = p 2 ( x ) - 4x 4 A 4 = F±(x)H 2 (x) 

or 

y 2 = {t 2 - Sl t + s 2 ) 2 - 4t 2 A 4 = (t 2 + at + b)(t + c) 2 , t = x 2 . 

There are three solutions for (si,s 2 ,b, c). The first is (^,0, ^ — 4A 4 ,0). To determine the 
classical limit we need to determine a. Similar to the discussion of 50(4) with one flavor we 
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need to consider various situations. The first is that there are D5-branes wrapping around the 
x = ±a, so a = —2a 2 , which gives the Special vacuum of 5*0(5) — > SO (3) x U(l) presented 
in the previous paragraph. The second is that all D5-branes are wrapped around the origin 
x = so we have 50(5) — > 50(5). Then we need a factor t 3 for the Chebyshev vacua where 
Nf + 2 = 3 and t° for Special vacua where 2Nf — N + 3 = 0. Since c = gives at least t 2 , we 
can only get two Chebyshev vacua by setting b = or a = ±4 A 2 , which match the counting 
(N — Nf — 2) = 5 — 1 — 2 = 2. The missing Special vacuum will be discussed immediately. 

The second and third solutions are (-a ± 2A 2 , (^±§^) 2 , t 2 ^) 2 , ^^-). For nonzero D5- 
branes wrapping around the x = ±a, again we have a = —2a 2 , so there are two vacua for 
50(5) — > U(2). Beside that, if we set a = 2A 2 for the second solution or a = — 2A 2 for the 
third solution, both will give the same curve having a complete square form, which is identical 
with the Special vacuum of 50(5) — > 50(5). 

We want to emphasize again that the calculations here are exactly the same calculations 
of 50(4) with Nf = 0. The only difference is how to explain and describe these results in the 
phase structure of 50(5) gauge group. 

5.3 S0(6) case 

For the 50(6) gauge group, the number of flavors Nf can be 0, 1,2, or 3. There exist the 
following four breaking patterns: 50(6) -> 50(2) x 17(2), 50(6) -> 50(4) x 17(1), 50(6) -> 
50(6), and 50(6) — > U(3). For the 50(2) factor, only the Special branch exists because for 
the Chebyshev branch there is a negative vacuum number N — Nf — 2 < 0, which does not 
exist when for 50(4) and 50(6) factors, both the Special branch and Chebyshev branch exist. 
The Nf — case has been discussed in section 4 under massive flavors, so we are left only with 
N f = 1,2,3. 

. IV) = 1 

1. Non-degenerated case 

The curve should be 



y 2 = p 2 ( x ) - 4r 6 A 6 = x 2 F 6 (x)H 2 (x). 



Using F 6 (x) = W'(x) 2 + bx 2 + d, H 2 {x) = x 2 + a, and P e (x) = x 2 (x 2 — ti)(x 2 — t 2 ), we found 
four solutions for (t = t\t<i,u = t\ + t 2 , a, b, d). The first is (0, a 2 , 0, 0, — 4A 6 ), which gives 
y 2 = x & (W'(x) 2 — 4A 6 ) and P&(x) = x 4 (x 2 — a 2 ). Classically it gives the Chebyshev branch of 
50(6) -> 50(4) x U(l) with the counting N-N f -2 = 4-1-2 = 1. The factor (x 2 f is exactly 
the one we need for 50(4) with Nf = 1 at the Chebyshev point by noticing (Nf + 2) = 3. 

The second solution has d — 0, so the curve is y 2 = ( w ^ ) 2 + b [x 2 (x 2 + a)] 2 . At the 
classical limit it becomes (0, a 2 , 0, 0, 0) which gives 50(6) — > 50(4) x U(l). Because of the 
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factor (x 2 ) 2 , it is the Special branch by noticing (2Nf — N + 4) = 2. The counting for the 
number of vacua is also consistent. 

The third and fourth solutions have d — and the classical limit becomes (a 4 , 2a 2 , —a 2 , 0, 0), 
which gives SO (6) — > 5*0(2) x U(2). Since d = 0, they are in fact the Special vacua with the 
counting two coming from the U (2) factor. 

To discuss the smooth transition, we can use the above solutions directly by taking various 
limits. However, since the solutions are so complex, it is not easy to see the results. However, 
it is easy to see that the first solution cannot be smoothly interpolated to the other three 
solutions. We will discuss the relationship of the other three solutions immediately by another 
method. 

2. Degenerated case and Special branch 

The curve should be factorized as 

y 2 = P 2 (x) - 4x 6 A 6 = F 4 (x)H 2 4 (x), 

or 

y 2 = P 2 (t) - 4t 3 A 6 = F 2 (t)H 2 2 (t), t = x 2 . (5.27) 
To solve the problem, we parameterize as follows: 

F 2 (t) = (t-a) 2 + b, H 2 (t) = t 2 + ct + d, P 3 (t) = t 3 - Sl t 2 + s 2 t - s 3 . 
There are four solutions. The first is given by 



-36 3b 2 3b (6 + 4A 4 ) 3 



4A 2 ' 2 16A 4 4' 6 64A 6 ' 

~ b A2 b a2 j (& + 4A 4 ) 2 

a = 4A^ + A ' C= 2A^ + A ' d =~^A^- 

There is only one sensible limit b/A 2 = (3 ^ where P 3 (t) -»• (t+ f ) 3 , so it gives SO(Q) -> U{3). 
The second and third solutions have 

si = 3(l±iV3) 



8A 2 ' 



which has also only one sensible limit 6/A 2 = (3 ^ and gives 5*0(6) — > £7(3). In fact, for all of 
these solutions we should consider two cases: one is b = and the other, 67^0. For 6 = 0, all 
three solutions give the same curve (5.27) as a square form with P 3 (t) = t 3 + A 3 . This vacuum 
is the Special vacuum of 50(6) — > 50(6). For b ^ 0, these three solutions give the three vacua 
of 50(6) — > U(3). Since 50(6) — > 50(6) has 6 = and a factor t at the curve, there is no 
smooth connection between 50(6) and U(3). 
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The fourth solution has 



4A 6 




8A 6 
b 2 



b ' 
4A 6 



S3 = 0, 



a 



16A 6 



b 



In fact, since d — and s 3 = 0, we can factorize out the t 2 factor and the curve (5.27) is reduced 
to the degenerated case of SO (5) without flavor. Thus, we immediately get the following results: 



Special vacua for 5*0(2) x U(2), one Special vacuum for 50(4) x U(l), and three Chebyshev 
vacua (N - N f - 2) = (6 - 1 - 2) = 3 for 50(6) -> 50(6). 
• N f = 2 

1. Non-degenerated case 

For simplicity we will use t = x 2 . The curve is given by 



Due to the factor t on the right hand side, we write Ps(t) = t(t 2 — s±t + s 2 ), Hi(t) = t + a, 
and F 3 (t) = t 3 + bt 2 + ct + d. There are three solutions for (s 1; s 2: a, c, d) as a function of 
6. The first is (-6/2, 0, 0, b 2 /A - 4A 4 , 0) which gives 50(6) -> 50(4) x 17(1). Since F 3 (t) = 
t(t 2 + bt + (b 2 /A — 4A 4 )), it is, in fact, the solution of Special or degenerated case. The other two 
solutions are (-6±2A 2 , ±(6^4A 2 ) 2 , |(6^4A 2 ), ±(6^4A 2 ) 2 , 0), which give 50(6) -> 50(2) x[/(2). 
Again, since d = they are the solutions of the Special or degenerated case. The counting 
number two comes from the U(2) factor. Finally, using the relationship of W'(x) and F 3 (t) we 
get b = —2a 2 . There is no smooth transition between these three solutions. It is noteworthy 
that there is no Chebyshev vacua of 50(6) — > 50(4) x U(l) because the counting leads to 
(N — Nf — 2) = (4 — 2 — 2) = 0. 

2. Degenerated case and Special branch 
The curve should be 



Let us consider the solutions case by case. First, if there is a factor t A , we have Ps(t) = t 2 (t — s), 
H2(x) = t 2 , and i^^) — (t — s) 2 — 4A 4 . There are two situations to be considered. If there 
are D5-branes wrapping around the x = ±a, we need to have F 2 (t) = (t — a) 2 + d, so s = a 2 , 
d = —4 A 4 and 50(6) — > 50(4) x U(l). It is the same solution presented in the previous 
paragraph. The power 4 = (2Nf — A^ + 4) of t is the expected character of 50(4) at the Special 
vacua. If there are no D5-branes wrapping around the x = ±a, we need a factor t from 7 2 (t) 
which is equal to set s 2 = 4A 4 . The curve has a factor t 5 which gives the Chebyshev vacua of 
50(6) — > 50(6). The counting two is also consistent with the fact that (JV — Nf — 2) = 2. 



(1) There is a smooth interpolation between 50(2) xU(2) and 50(4) x U(l); (2) There are two 



y 2 = Pi(t)-4AH* = tF 3 (t)H 2 1 (t). 



y 2 = P 2 (t)-4AH* = F 2 (t)H 2 2 (t). 



(5.28) 
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If there is only one t 2 factor in (5.28), we write P 3 (t) = tP 2 (t) and H 2 (t) = t(t + a). 
Canceling out the factor t 2 , we get P 2 (t) — 4A 4 t 2 = F 2 (t)(t + a) 2 . It can be solved if we require 
P 2 (t) - 2i]A 2 t — (t + a) 2 , so F 2 (t) — (t + a) 2 + 4i]A 2 t. If there are D5-branes wrapping around 
x = ±a, we get a = —a 2 and two Special vacua for 5*0(6) — > SO (2) x U(2), as found in the 
previous paragraph. If there is no D5-brane wrapping around x = ±a, there are two choices: 
either F 2 (t) has a factor t so that a = or F 2 (t) is a square. The former where a = gives a 
factor t 5 , which is the one discussed before. The latter is given by P 2 (t) = t 2 + A 4 , which gives 
the one vacuum of 5*0(6) — ► 50(6) at the Special branch. 

Finally we consider the case in which there is no t 2 factor at all, so all D5-branes must wrap 
around x = ±a. To achieve this, we must have 

P 3 (t) - 2A 2 t 2 = (t + a) 2 {t + b), P 3 (t) + 2A 2 t 2 = (t + cf{t + d) . 

This problem is the same problem of U (3) with Nf = 4 which has been solved in [20] 

,2 2siA 2 , ^ a2 n 2s x A 2 
a = si + s 2 , C — S1—S2, b = si - 2A - 2s 2 H , d = s 1 + 2A + 2s 2 H , 

s 2 s 2 

with s\(s 2 + A 2 ) = s 2 A 2 . Furthermore, using (t + b)(t + d) = (t — a 2 ) 2 + e we get 

1 2A 2 ^ 



bd = -2a z -> si ^1 + —J = -a 2 - 

For fixed a, there are three solutions which have the classical limits (si, s 2 ) — > (—a 2 , 0). These 
three vacua give 50(6) — > C/(3). 

It is worth remarking about the different behaviors. We observe a smooth interpolation 
between 50(4) x U(l) and 50(2) x U{2) for iV> = 1 at Special branch, but not for N f = 2. 
The reason for that comes from the power of t. The Special branch of 50(2) always has a factor 
t 2 in the curve. For 50(4), it is (2Nf — N + 4) = 2Nf, so only for the case with Nf = 1 can 
we have a factor t 2 . This also explains why there is no smooth interpolation of 50(4) x U(l) 
or 50(2) x U{2) to U(3) because there is no t factor for the curve of U(3). 

. JV> = 3 

1. Non-degenerated case 

Setting t = x 2 , the curve of non-degenerated case should be 

y 2 = P 2 (t)-4A 2 t 5 = tF 3 (t)H 2 (t). 

Writing P 3 (t) = t(t 2 - stf + s 2 ), H^t) = t + a, and F 3 (t) = t 3 + bt 2 + ct + d, we found three 
solutions. The first one has (si, s 2 , a, c, d) = ( ~ b ~ 4A2 , 0, 0, ^ 6+4 4 A ^ ). It gives one vacuum of 
50(6) — >■ 50(4) xU(l). However, since F 3 (t) = t(t 2 + bt + c), it is in fact at the Special branch. 
It is also noteworthy that the power of t is 4, which is not equal to (2Nf — N + 4) =6. The 
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reason is that SO [A) with three flavors is not asymptotically free (recalling that Nf < N — 2); 
thus, the power stops at N. In fact, we have a similar result that there is a factor t 2 for 5*0(2) 
with flavors. The other two solutions have d — and the classical limits (—6, b 2 /A, 6/2, 6 2 /4, 0), 
so give 50(6) — > 50(2) x U(2). Again, since d — 0, these two vacua are in fact at the Special 
branch. The counting number two comes from the U (2) factor. Finally, we can fix as follows: 
b = — 2a 2 — 4A 2 . The reason why we did not find the Chebyshev branch in the non-degenerated 
case is that the number of counting becomes (N — Nf — 2) < for both the 50(2) and 50(4) 
gauge groups. 

2. Degenerated case and Special branch 
Let us now consider the Special branch or the degenerated case with the following curve 

y 2 = P 2 (t)-4A 2 t 5 = F 2 (t)H 2 (t). 

We will discuss them case by case. First, if there is a t 4 factor in the curve, we must take 
P 3 (t) = t 2 {t - s), H 2 (t) = t 2 and F 2 (t) = (t - s) 2 - AA 2 t. There are three possible values for s. 
If s = —A 2 , the curve is a complete square and gives the Special vacuum of 50(6) — > 50(6). 
Another case is s = a 2 which gives one Special vacuum of 50(6) — > 50(4) x U(l). As we have 
remarked before, since Nf = 3 > (N — 2) = 1, the power of t is N — 4. Finally, if F 2 (t) has just 
one factor t, or if s — 0, it gives the one vacuum of the Chebyshev branch 50(6) — > 50(6). 

Second, if there is only t 2 factor in the curve, we have Ps(t) = t(t 2 — sit+s 2 ), H 2 (t) = t(x+c). 
Cancelling the t 2 at both sides, we get the reduced curve 

P 2 (t) -4A¥ = F 2 (t)H 2 (t). 

This is exactly the same curve (5.26) and we get two Special vacua of 50(6) — > 50(2) x U(2) 
by keeping only those solutions with a factor t 2 . 

The last case is when there is no factor t in the curve and all D5-branes wrap around the root 
x = ±a, so the gauge group is broken to 50(6) — > U(3). Using this fact, we know immediately 
F 2 (t) = (t - a 2 ) 2 + d - 4A 2 t. There are three solutions 24 . 

The phase structure of 50(6) with massless flavors has been summarized in Table 7. 

5.4 S0(7) case 
The curve is 

y 2 = P 2 (t)-AA l0 - 2N ft 1+N f. 

24 It is very complicated to solve the factorization. However, notice that at the limit A — > 0, there are 
a — > a 2 , c — ► —2a 2 , d — > a 4 and 6^0. If we set a = 1 and A = 0.001, numerically it is easy to see only three 
solutions that satisfy this limit 
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3 





Table 7: Summary of the phase structures of 50(6) gauge group with massless flavors. It is 
worth comparing with the the phase structures of 50(6) with massive flavors. Here again we 
use capital letters in the last column to indicate the phases having smooth transition. 
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For Nf = 0, it has been discussed in the section 5 of massive flavors. For Nf > 1, the curve 
can be written as 



y 2 = p|(t)-4A 10 - 2JV /t 2 t JV /- 1 

which is the same curve of 5*0(6) with (Nf — 1) flavors. Therfore, all the results of 5*0(6) can 
be applied here. 

5.5 The smooth interpolation 

Now let us discuss in this subsection the general picture of smooth interpolations. The 
classical limit has three types: SO(N c ) -> SO(N ) x U(N X ), SO(N c ) -> U([N c /2\), and 
SO(N c ) — > SO(N c ). The first two have U(l) left at IR while for the last one, there is no 
U(l) left at all. Because of this fact, only the first two types can have smooth interpolations 
by realizing the presence of the number of U(l) left. 

Now for the type SO(N c ) — > U([N c /2]), the curve has no factor t in the factorization form, 
while for the type SO(N c ) — > SO (No) x U(Ni), depending on the Chebyshev branch or Special 
branch, it has a different power of t. 

For the Chebyshev branch, the power of t = x 2 is always an odd number and mainly the 
function of Nf. Thus, the type SO(N c ) — > SO (No) x U(N±) cannot be smoothly connected to 
the type SO(N c ) — > U([N c /2}), but can have smooth interpolations inside itself. For example, 
if SO(N c ) -> S6(N ) x U(N^ and SO(N c ) -> SO(Mo) x U(M x ) with N f < N - 2 and 
iVj < M — 2 (so the Chebyshev branch exists for both SO(N ) and 50 (M )), the curves will 
have the same power of t and can give the above two different classical limits. To demonstrate 
the analysis, let us consider the 50(8) with one flavor. The curve is 

y 2 = P£(x) - Ax 6 A 10 = x 2 F 6 (x)Hl(x). 

The required x 6 factor forces us that Ps(x) = x 4 (a; 4 — s\x 2 + s 2 ), H 4 (x) = x 2 (x 2 + a). Writing 
F 6 (x) = x 6 + for 4 + cx 2 + d, we have the following solutions 

iA 5 , 2 iA 5 A 10 , A 10 



>A 5 2 „ iA 5 



Sl = _2a±— — , s 2 = a =f3 



a s/2' ^ - -r " a i/ 2 - 

We can take the following limits: (1) A -> 0, a fixed, we get 50(8) -> 50(4) x U(2); (2) A -> 0, 
a — > 0, but ^ — > /3 7^ 0, we get 50(8) — > 50(6) x U(l). We have not shown that it is always 
true that SO(N Q ) x U(N 1 ) is connected to 5O(M ) x f/(Mx) under the above conditions, but 
we expect it to be true. 
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Now let us consider the Special branch in which the power of t = x 2 is always an even 
number. 1) First, if the power is zero, it is possible to connect SO(N c ) — ► 5O(iV ) x U(N 1 ) to 
SO(N c ) -> U([N c /2\). Recalling that the power is given by (2N f -N+4)/(2N f -N+3) for N is 
even/odd and N f < N-2 and N/(N-1) for N f > N-2, we need (2N f -N+4)/(2N f -N+3) = 
or N f = (N - 4)/2 for N even and N f = (N - 3)/2 for N odd. This has been observed for 
SO (5) with Nf — where there is, in fact, a smooth interpolation in 5*0(3) x [/(l) <-> t^(2). 

2) Second, if the power is not zero, we can only expect a smooth interpolation inside the 
first type. Since the power is a function of both Nf and N, the symmetry breaking pattern 

5O(iV ) x U(Ni) 

is connected to the following symmetry breaking pattern 

SO(M ) x [/(Mx) 

(assuming M < N and both are even numbers) only if Nf > M — 2, Nf < N — 2 and 
(2N f - iVo + 4) = M . One such example is S0(6) with one flavor where S0{2) x C/(2) is 
smoothly connected to 5*0(4) x U(l). In fact, case 1) where the power of t is zero and case 2) 
where the power of t is not zero are related to each other by the addition map. 

It is noteworthy that for the smooth transition in the Special branch, the condition (2Nf — 
Nq + 4) = M is exactly the condition of the Seiberg dual pair between SO(N ) and 50 (M ). 
Therefore, the smooth transition can be connected through the Seiberg duality. However, for a 
smooth transition in the Chebyshev branch, no such relationship exists and a smooth transition 
cannot be understood from the Seiberg duality. 
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Appendix A Strong gauge coupling approach: superpotential and a 
generalized Konishi anomaly equation for SO(2N c ) case 

In [19] the strong gauge coupling approach for SO/USp pure gauge theories were studied 
generally. That analysis was an extension of [8] to allow a more general superpotential in 
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which the degree can be arbitrary without any restrictions. We extend the analysis discussed 
in [19] to SO(2N c ) gauge theory with N f flavors by using the method of [72, 79, 80]. Let us 
consider superpotential regarded as a small perturbation of an Af = 2 SO(2N c ) gauge theory 
[65, 66, 67, 68, 69, 70, 71, 72, 8, 11] 

W(*) = E f 1 ^ 28 = E 92 S u 2s , u 2s = -Tr$ 2s (A.l) 

s=l ZS s=l ZS 

where $ is an adjoint scalar chiral superfield and we denote its eigenvalues by ±0/(7 = 
1, 2, • • • , N c ). The degree of the superpotential W(§) is 2(k + 1). Since $ is an antisymmetric 
matrix we can transform to the following simple form, 

$ = (_io) ® dia g(^i'^2,-"-,^ c -r,0,0,---,0). (A.2) 

When we replace Tr& 2j with (Tr& 2j ), the superpotential becomes an effective superpotential. 
We introduce a classical 2N C x 2N C matrix $ d such that (Tr$ 2 -?) = Tr$ 2 / for j = 1,2,---, N c . 
Additionally, u 2j = f-Tr$ 2 / are independent. However, for 2j > 2N C , both Tr& 2j and (Tr$ 2i ) 
can be written as the u 2 j of 2j < 2N C . Classical vacua can be obtained by putting all the 
eigenvalues of <£> and <3> ci equal to the roots of W'(z) = Ylt=\ 92rZ 2s ~ l ■ We will take the degree of 
superpotential to be 2(k + l) < 2N C first in which the u 2j are independent and (Tr& 2j ) = Tr$ 2 /. 
Then we will take the degree of superpotential to be arbitrary. Until now we have reviewed the 
discussion given in [19]. Next we will study some restriction that is an important idea in the 
discussion below. 

Let us consider the special (A^ c — r) dimensional submanifold of the Coulomb branch 
where some of the branching points of the moduli space collide. The index r runs from 
to min(A r c , Nf/2) [68, 80] and classifies the branch. On the r-th branch, the effective theory 
becomes SO(2r) x U(l) Nc ~ r with Nf massless flavors. Thus, after turning on the tree level 
superpotential since there exist U{l) n gauge groups with 2n < 2k, the remaining (N c — r — n) 
U(l) factors are confined, which lead to (N c — r — n) massless monopoles or dyons. This occurs 
only at points where at W = the monopoles are massless on some particular submanifold 
< u 2r >. This can be done by including the (N c — r — n) monopole hypermultiplets in the 
superpotential. Then the exact effective superpotential by adding (A.l), near a point with 
(N c — r — n) massless monopoles, is given by 

N c -n-r k+1 

W eff = V2 E M i(u 2s )qiqi + E 92 S u 2s . 

1=1 s=l 

By varying this with respect to u 2s , we get an equation of motion similar to pure Yang-Mills 
theory except that the extra terms on the left hand side since the u 2s with 2s > 2(N C — r) 
are dependent on u 2s with 2s < 2(N C — r). There exist (N c — n — r) equations for the (k + 1) 
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parameters g 2s . Here qi and q\ are the monopole fields and Mi(u 2s ) is the mass of Z-th monopole 
as a function of the u 2s . The variation of W e ff with respect to qi and q% vanishes. However, 
variation of W e ff with respect to u 2s does not lead to the vanishing of qi% and the mass of 
monopoles should vanish for / = 1, 2, • • • , (N c — r — n) in a supersymmetric vacuum. Therefore, 
the superpotential in this supersymmetric vacuum becomes W exact = Y%=1 9?s < u 2s >■ The 
masses Mj are equal to the periods of some meromorphic one-form over some cycles of the 
N = 2 hyperelliptic curve. 

It is useful and convenient to consider a singular point in the moduli space where (N c — r — n) 
monopoles are massless. Then the Af = 2 curve of genus (2N C — 2r — 1) degenerates to a curve 
of genus 2n [8, 11] and it is given by 25 

= p2^_ r){x) _^ { N c -r)-2 { N } -2r + 2) x 2 { N s -2r + 2) 

= ^ 2 ^2iV c -2r-2n-2(^) i71 2(2n+l)(a:), (A.3) 

where 

N c -n-r-l 2n+l 

H 2Nc - 2r - 2n - 2 (x) = ]J (x 2 -p 2 ), F 2{2n+1) (x) = (x 2 - q 2 ). 

i=i i=i 

Here H 2Nc _ 2n _ 2r - 2 (x) is a polynomial in x of degree (2N C — 2r — 2n — 2) that gives (2N C — 
2r — 2n — 2) double roots and i 7 2(2n+i)(^) is a polynomial in x of degree (An + 2) that is related 
to the deformed superpotential. Both functions are even functions in x. That is, a function 
of x 2 which is peculiar to the gauge group SO(2N c ). Since we concentrate on the (iV c — r) 
dimensional subspace of Coulomb phase the characteristic function P 2 (N c -r){ x ) is defined by 
2(N C - r) x 2(N C - r) matrix $ d as follows: 

N c -r 

P2 {Nc -r)(x) = det(:r - $ ci ) = XKx 2 -4>)). (A.4) 

The degeneracy of the above curve can be checked by computing both y 2 and at the point 
x = ±pi and x = obtaining a zero. The factorization condition (A.3) can be described and 
encoded by Lagrange multipliers [18]. When the degree {2k + 1) of W'(x) is equal to (2n + 1), 
we obtain the result (A. 6). 

• Field theory analysis for superpotential 

In [20] it was noted that when the degree of W'{x), n, was greater than or equal to (2N C —Nf) 

the structure of the matrix model curve of U (N c ) gauge theories was changed by the effect of 

the flavors. Thus, we first extended the discussion in [20] to SO(2N c ) gauge theory and then 

25 As we classified in previous section, there exist various curves (3.4), (3.5) and (3.6) depending on the number 
of single and double roots. It is straightforward to proceed to other cases so we restrict ourselves to here the 
particular case (A.3). 
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later to the more general cases. The massless monopole constraint for SO(2N c ) gauge theory 
with Nf flavors is described as follows: 



V 2 = P| w Jx)-4A 4 ^- 4 - 2 ^A(x)=x 2 ^_ 2n _ 2 ( a ;)F 2{2n+1) (a;), 

l-l 2 N f 

= x 2 n(^ 2 -P-) F 2{2n+1) (x), A(x) = x 4 H(x 2 - m)) (A.5) 

i=l j=l 

where we used / as the number of massless monopoles. From the equation (A.5) we can find 
double roots at x = 0, ±pi i = 1, — 1. We do not need take into account all points. 
Since P 2 n c (x) and A(x) = x 4 detjv / (x + m) are an even function of x, we have only to consider 
x = 0,x = +pi. For simplicity of equations we introduce pi = 0; thus, index i runs from 1 
to I. With these conventions we have an effective superpotential with / massless monopole 
constraints (A.5), 

W low = Y j g2tU2t + jz\u(P2N c {Pr)-^ Nc ^ Ni 
t=l i=l L ^ 

+ Qi^~ (P2NM) - 26 l A 2 ^- 2 -^ V / Ife)) 



where Lj and Qi are Lagrange multipliers and t; L = ±1. From the equation of motion for p^ and 
Qi we obtain the following relations, 

Q l = 0, (p2nM) - 2eA 2Nc - 2 - N{ \[M&)) = 0. 

The variation of W\ ow with respect to u 2 t leads to 

9* + ( P ^ c fe) - 2e l A 2 ^- 2 "^ V / Ife)) = 0. 

Since A(pi) is independent of u 2 t the third term vanishes. By using P2N c iPi) = J2f=o s 2jP 2 Nc ~ 2 ^ 
26 we obtain 

l Q i n c 

92t = J2 L i* P ™ c (Pi) =J2J2 L iP 2 i NC ~ 23s 2j-2t- 

i=l OU 2t i=l j= 

With this relation as in [12] we can obtain the following relation for W'(x). 

xPom(x) 



W'(x) = Y,^f^^- x ~ 1 T,^L t A 2N ^- N f^A( Pl )+0(x- 3 ). 

i=l X Pi i=l 



26 The polynomial P 2 N c (Pi) is given by J2f=o s 2jX 2Nc 2j = ]l/=i( x2 — <Aj) where s 2j and u 2j are related each 
other by so-called Newton's formula s 2 j + X)i=i s 2j-2iU 2 i = where j = 1,2, ■ ■ ■ , N c with sq = 1. From this 



recurrence relation we obtain = —s 2 j-2t f°r j > t. 
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Defining a new polynomial -82(2-1) of order 2(1 — 1) as in [2], 



E 



Li _ B 2 (i-i){x) 



{x 2 -p 2 x 2 H 2 (i~i){xY 



the W'(x) is rewritten as 



W\x) = xP ™£ ) ^ l 7 1 l ix) +Q(x- 1 ). 
x 2 H 2{ i^ 1) (x) 

Now we can compare them with both sides, in particular, for the power behavior of x. It is 
easy to see that the left hand side behaves like (2n + 1) while the right hand side behaves like 
2N C — 1 — 2(1 — 1) except the factor B 2 n_u(x) and therefore the condition I = (N c — n) will 
give rise to the consistency and the polynomial B 2 n_i\ (x) becomes constant. By using this 
relation and substituting the characteristic polynomial P 2 n c (x) from the monopole constraints 
(A. 5), we can obtain the following relation together with the replacement of the polynomial 
H 2 (i-i)(x), 

A T 2 \4N c -2Nf-4 yi N f I „2 _ m 2\ i 

W*) + U N c - n J, 2v ^ ~ = K +1 (*) 2 + 0(x*j) . (A.6) 

llj=l \ x ~Pj) i/2n+2 

Thus, if n > —Nf + 2N C — 3, the effect of flavor changes the geometry. This is a new feature 
compared with pure gauge theory with no flavors. This is due to the fact that the flavor- 
dependent part, A(x) provides a new contribution. 

The formula (A.6) makes the assumption that every root of W'(x) has some D5-branes 
wrapping around it. If the assumption is not true, for example, for the breaking pattern 
SO(N c ) -> n"=i U(Nj), instead of SO(N c ) -> SO(N ) x U]=i U(Nj), the formula (A.6) should 
be modified. The general form can be obtained as follows. Assuming a curve is factorized as 

y 2 = P 2 M (x) - A(x) = f 2s (x)H 2 M _ s (x) (A.7) 

where A(x) counts the contribution of flavors (massive or massless) and Hm~ s (x) includes all 
the double roots (for example, the x 2 in [A. 3] should take into account H M _ s [x]), we will have 

/*(*) + n A{X L 2 = K(x) 2 + 0(x s - 1 ) (A.8) 

where W' s (x) = Y\ S j =1 (x — Xj) and the A/s are the roots of W'(x) with wrapping D5-branes. For 
example, for the case considered in (A.6), Xj = 0, ±aj so W' s (x) = W^n+i^)- ^ among these 

W (x) 

roots there is only the origin without wrapping D5-branes, W' s 



F4n(x) + TrJVc-n/o J2\ 2 ~ „2 



Uf=r(x 2 -P 2 ) 2 9ln + 2 
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X 



2 

i +0(x l 



W 2n+ i(x) \ ;>„-> . 



(A.9) 



As soon as n > —Nf + 2N C — 3, the effect of flavor also changes the geometry. 
• Superpotential of degree 2(k + 1) less than 2N 

Now we generalize (A. 6) to 2n < 2k by introducing the constraints (A. 3). We follow the 
basic idea of [18, 19] and repeat the derivations of (A. 6) and generalize to the arbitrary degree 
of superpotential later. At first, in the range 2n + 2 < 2k + 2 < 2N, let us consider the 
superpotential for SO(2N c ) theory with massless Nf flavors in the r-th branch for simplicity 
(the massive flavors can be obtained by including the mass dependent parts appropriately. For 
massless case, as we remarked before the r is fixed to be some number for Chebyshev or Special 
branches.) under these constraints (A. 3) by starting from the curve in pure case and adding 
the flavor dependent parts, 



fc+l 2N c -2r-2n-2 

W eff = Y,92sU 2s + Yl 
s=l i=0 



{x-piY 



X-pi) 

(A.10) 



where Lj and Bi are Lagrange multipliers imposing the constraints and = ±1. What we have 
done here is to include the r-dependence and Nf dependence coming from the generalization of 
pure gauge theory [19]. For an equal massive case, we simply replace x Ns ~ 2r+2 in the numerator 
with (x 2 - m 2 ) N f/ 2 ~ r x 2 . The cont our integration encloses all pi's and the factor l/2m is 
absorbed in the symbol of § for simplicity. The pi's where % — 0, 1, 2, • • • , (2N C — 2r — 2n — 2) 
are the locations of the double roots of the curve y 2 = Pi(N c -r)( x ) — 4A 4Arc ~ 2Ar /~ 4 a? 2Ar - f ~ 4r+4 
according to the constraints (A. 3) where the function H 2 N c -2r-2n-2(x) has a factor (x 2 — p 2 ). 
The function P 2 (N c -r)( x ) depends on u 2s . Note that the massless monopole points appear 
in pair (pi, —pi) where i = 1, 2, • • • , (N c — r — n — 1). Therefore, we denote half of the pis 
by VN a -r-n-i+i = — Pi, i = 1, 2, • • • , (N c — r — ii — 1). Moreover, we define p = and note 
that the summation index % in the above starts from % — 0. Since B 2 (]s[ c -r){ x ) is an even 
function in x and the property of if the following constraints are satisfied at x = pi where 
i = l,2,...,(7V c -r-n-l), 

(P 2 (^^)(^)-2^- 2r+2 QA 2JVc ^^ 2 )U= Pl = 0, 
^(p 2( ^ c _ r) ( a ;)-2^"^ 2 e l A 2 ^^- 2 )U =Pi = 0, (A.ll) 

then they also automatically are satisfied at x — —pi- Then the numbers of constraints that we 
should consider are (A^ c — r — n — 1). Thus, we denote the half of the Lagrange multipliers by 
LN c -r-n-i+i = and -Btv c -7— n-i+i = where % — 1, 2, • • • , (A^ c — r — n — 1). Due to the fact 
that the second derivative of y with respect to x at x = pi does not vanish, there are no higher 
order terms such as (x — Pi)~ a ,a = 3,4, 5, • • • in the effective superpotential (A. 10). In other 
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words, for given constraints (A. 3), there exist only a Lagrange multiplier Lj with (x — Pi)^ 1 
and a Lagrange multiplier Bi with (x — Pi)~ 2 - 

The variation of W e ff with respect to B { leads to an integral expression which is the coeffi- 
cient function of Bi in (A. 10) and by the formula of contour integral, it is evaluated at x = Pi 
and we have changed the derivative of P2(N c -r)( x ) with respect to x into P 2 (N c -r)( x ) an d the 
trace of some quantity 

[x-piY 



= (P 2{Nc - r) (x) - 2x N f-^e l A 2N ^- 2 )' 



X=Pi 



N f -2r + 2 i 



Ob 



= (P 2 (iV c -r)(z) X, X 2 _ g ^2(AT c - r )(^) 

1 ^ + 2^ 



|X=Pi 



= P 2 (^ c -r)(a;) Tr 



X-$ cl X 

where we used the equation of motion for L { when we replace 2x N f~ 2r+2 eiA 2Nc ~ N f~ 2 with 
P 2 (N c -r)( x ) at x — Pi. 27 The equation of motion for Bi can be summarized as 28 

(Tr * - / + 2 ) \ x=Pi = 0, P 2 {N c -r){x = Po = 0) = 0, P 2(i v c -r)(a: = Pi) ^ 
where the characteristic function by solving the first order differential equation is given by 

/ oo 



P2(N c -r)(x) 



X 



«-"exp(-g^| 



(A.12) 



and the polynomial part of a Laurent series inside the bracket is denoted here by +. By 
putting the negative power of x to zero, the u 2s with 2s > 2(N C — r) can be obtained the u 2s 
with 2s < 2{N c -r). 

Next we consider the variation of W e ff with respect to pi, 

f Pw -r)(x) - 2x N f~ 2r+2 e t A 2Nc ~ Nf - 2 , 

= 2B,. & 2(Nc r)K ' — - — dx 

J (x — Pj) 6 



27 The last equality comes from the following relation, together with (A. 2), Tr ^ 1 ^ = J2T=o x k lr ^ r ^cl = 
ESo^" (2i+1) ^ci = YZ^^ 2 Y {l+1) T.t^m l = (Ef=r^) + f where $ c , is antisymmetric 
matrix, and the odd power terms are vanishing. 

28 Through the definition (A. 4), the derivative P 2 (iv c -r)( a; ) with respect to x is given by P2(N c -r)( x ) = 

(j\f c - r (x 2 - cbj))' = 2xj:^- r l\^7( x2 - </>?) - P2(N c -r){x) Ej=7 Taking into account the result 
of trace given in previous footnote 27, we can rewrite this result as p 2(JVc ~'' ) ^l = Tr— ^ ^r- Therefore, we 

have p^ c I^(lj = T l Zo x ' {2l+1)Tl ^li = J2Zo x ' i2l+1)2iu ^- lt is eas y to check that for °q ual massive flavors, 
the corresponding relation leads to (Tr^-±- - J N ^ - \ x=Pi = 0, P 2 (N c -r)( x = Po = 0) = 0, P 2 (N c -r)( x = 
Pi) + 0. 
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where there is no L; L term because we have used the equation of motion for In general since 
this integral does not vanish, we should have Bj = because the curve does not have more than 
cubic roots due to the fact that y 2 contains a polynomial H 2Nc _ 2r _ 2n _ 2 (x) as we have discussed 
before. 

Let us consider variation of W e ff with respect to u 2s , by using the relation dI>2 ^- r )^ — 



p. 



2(N c -r)(x) 



, which can be checked from the definition (A. 12) and remembering the fact that 



the function P 2 (A r c -r)(^) depends on w 2 s, 

P2(N c -r)(x) 



2N c -2r~2n-2 
= 92s - E 



i=Q 



x 2s 



— dx, 

X-Pi 



where we used -Bj = at the level of equation of motion. Multiplying this with z 2s 1 and 
summing over s where z is inside the contour of integration, we can obtain the first derivative 
W'(z), 

k+1 2N c -2r-2n-2 k+1 p ( \ r 

w(z) = s>^-> = £ fT,^- l 2(K -2 ,{> ,,1 (A.i3) 

s=l i=0 s=l x \ x Vi) 

Let us introduce a new polynomial Q(x) defined as 

2N c -2r-2n-2 r N c -r-n-l 9 2r f)/^ 

E T^ = L 0+ E J T~ b 2 = rr ) 1 s (A-14) 

where we used the fact that Ljv c -,-- n -i+i = Li and pN c ~ r -n-i+i = ~Pi for i = 1, 2, • • • , (A^ — 
r — n — 1). By using this new function we can rewrite (A.13) as 

Since W'(z) is a polynomial of degree (2k + 1), we found the order of Q(x) as (2k — 2n), so we 
denote it by Q2k-2n(x)- Thus, we have found the order of polynomial Q(x) and therefore the 
order of integrand in (A.13) is like 0(x 2k ~ 2s+1 ). Thus, if s > k + 1 it does not contribute to 
the integral because the power of x in this region implies that the Laurent expansion around 
the origin vanishes. We can replace the upper value of summation, k + 1, with the infinity and 
by computing the infinite sum over s we get 

^(^)=/E^ g2 " 2w(3:)P2(jVC "; )( f ) ^=^ Q^n(x)P2 {Nc - r) (x) 

J ^ X 2s xH 2Nc -2r-2n-2(x) J x(x 2 - Z 2 )H 2Nc ^ 2r _ 2r ^ 2 (x) 

From (A. 3) we have the relation, 



P2 {Nc -r)(x) = x^F 2i2n+1) (x)H 2Nc _ 2r _ 2n _ 2 (x) + 0(x- 2N < +2N f- 2r+4 ). (A.16) 

Let us emphasize the presence of the second term which will play the role and cannot be ignored 
because the power of x can be greater than or equal to —1 depending on the number of flavor 
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Nf which can vary in various regions. For pure cases, there were no contributions from the 
second terms. Therefore, we have 

W\z) = f z^f^dx + f O (x~ m ^ +2N f +2k+3 ) dx, (A.17) 
where we have defined y m as follows, 

Vm( X ) = F 2{2n+l)(x)Ql k _ 2n (x) (A.18) 

corresponding to the matrix model curve. In (A.17), the second term does contribute if the 
condition — AN C + 2Nj + 2k + 3 > — 1 holds. Then we get an expected and generalized result: 

2 / \ t? nrf ,,_{ W 2 ' k+ 2 {x) + 0(x 2k ) k>2N c -N f -2 

vj?) - F 2(2n+1) (x)Q 2k _ 2n (x) - j w , k+i 2 (x)+a ^ Nc -2 Nf -^ k<2Nc _ Nf _ 2 

= W 2 ' k+1 2 (x) + f 2M (x), 2M = mzx(2k,4N c -2N f -4) (A.19) 

where both F 2 (2 n +i)(^) and Q2k-2n(x) are functions of x 2 , then /W^) is also a function of x 2 . 
We put the subscript m in the y m in order to emphasize the fact that this corresponds to the 
matrix model curve. When 2k = 2n, we reproduce (A. 6) with Q = Q2n+2- The second term 
on the left hand side of (A.6) behaves like as x ^N s -i{N c -n-i) = x 2N f -4N c +4n+e _ Depending on 
whether the power of this is greater than or equal to 2n = 2k, the role of flavor is effective or 
not. When n = k > —Nf + 2N c — 3, the flavor dependent part will contribute to the W'(x). The 
above relation determines a polynomial F 2 ^ 2n+ i^(x) in terms of (2n + 1) unknown parameters 
by assuming the leading coefficient of W(x) to be normalized by I. These parameters can be 
obtained from both P2(N c -r)(x) and H 2Nc ^ 2r _ 2n _ 2 (x) through the factorization condition (A. 3) 

29 

When k is arbitrary large, we refer to the version one in the hep-th archive for details. 
• A generalized Konishi anomaly 

Now we are ready to study the derivation of the generalized Konishi anomaly equation 

based on the results of previous section. As in [18], we will restrict ourselves to the case with 

29 One can proceed the case where there are no D5-branes wrapping around the origin similarly. Our 
corresponding monopole constraint equation (A. 16) becomes P2(N a -r)( x ) = \J F4 n (x)H2N c -2r-2n{x) + 
0{x~ 2Nc+2h Tj : ~ 2r + 4 ) . Using this relation, it is straightforward to arrive the following matrix model curve in 
this case: 

f^k±l^l) 2 + (x 2k - 2 ) k>2N c -N f -2 
y 2 m {x) = F in {x)Q 2 k ^ n {x) = { ) x K 

+0 k<2 N c -N f -2 



X 



f 2M {x), 2M = max(2fc-2,4A^ c -2iV / -6). 
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(TrW(<E>)) = TrW'($ d ) and assume that the degree of superpotential (2k + 2) is less than 27V C . 
By substituting (A. 14) into (A. 15) we can write the derivative of superpotential W'((f)j) 



2N c -2r-2n-2 

(x 2 - <pj) (x - Pi) 



i=0 



where we varied W(<f>i) with respect to 0/ rather than -u 2s and used the result of Bi = 0. Note 
that the characteristic function is given by P 2 ( Nc - r )(x) = Iljiiv^ 2 — 4>j). Using the above 
expression, one obtains the following relation 30 , 

W'(® ,) Nc ~ r 1 



Z-^cl S (z 2 -(j)j) 



s (^ 2 - h J (* 2 - (* - ft) 

2JVC "E 2n " 2 7 ft^^ , f^Tr^— - xTr^— ) d*. (A.20) 
i=o ( x - 2 )( x - ft) ^ ^ - $ d a: - 

As in the case of [18], we can rewrite outside contour integral in terms of two parts as follows: 



/ (A.21) 

JC,+C-, 



where C z and C_ z are the small contour around z and — z respectively. Thus, the first term in 
(A.20) (corresponding to the second term in [B.3] of [18]) can be written as, by exploiting the 
relation (A. 14) to write in terms of Q2k-2n(x) and H 2 N c -2r-2n-2(x) 

1 \ I zQ 2k -2n(x)P2(N c -r)(x) 



Z-®cl) Jz out xH 2 N c -2r-2n-2(x)(x 2 - Z 2 ) 

Let us emphasize that in this case also we cannot drop the terms of order Q( x - 2N c+ 2N f- 2r +^ 
in the characteristic function P 2 (N c -r)(x)- In order to take into account of this, we have to use 
the above change of integration, then the first term of (A.20) is given by 

zQ2k-2n(x)P2(N c -r)(x) _ f zQ 2k ~2n(x)P2(N c -r)(x) ^\ 




xH 2 N c -2r-2n-2(x)(x 2 - Z 2 ) Jc z +C- Z xH 2 N c -2r~2n-2(x)(x 2 - Z 2 ) 

ym(z)P 2 (N c -r)(z 



P i(N c -r)( Z ) - 42 2 ^- 4r + 4 A 4 ^- 2 ^- 4 



(A.22) 



30 When we change the summation index from k to i, the only odd terms appear because effectively the prod- 
uct of $ c/ and W'($ c j) does contribute only under that condition Tr ^§^ = TrJ^kLo z^' 1 § k cl W c i) = 
E l =o z~ (2i+1)_1 2 Ef=r <t>Y +1 W{<t>i) = 2 Ef=r <l>iW{<f>i) , z2 ^ r The even terms do not contribute. Here z 
is outside the contour of integration. We recognize that the following factor can be written as, by simple manip- 
ulation between the property of the trace we have seen before, ( z 2_ l p^ x 2_^i) = ( x i- z -i) i^^ 1 z \ \ ~ x ^ r x ; ) ■ 
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where the first term was obtained by the method done previously and the second term was 
calculated at the poles and we used 

J P i(N c -r)( Z ) - 4A^- 2N f-*Z 2N f- 4r + 4 

H 2Nc _ 2r _ 2n - 2 (z) = , Vm(z) = F 2(2n+\) (z)Ql k - 2n ( Z ) ■ 

Z yj F 2(2n+1) (z) 

The crucial difference between U(N C ) case and SO(2N c ) case comes from the second term 
of (A. 20) (corresponding to the first term in (B.3) of [18]), which vanishes in U(N C ) case. 31 
Now we use the result of the equation of motion for Bi (the equation just above (A. 12)) in 
order to change the trace part and arrive at the final contribution of the second term (A. 20) as 
follows: 

2JV c -2r-2ri-2 r p / _ \ 

t^o (Pi ~ z ) 

2 ^-g 2 - 2 f (N f + 2)P 2 ^- r) (x) U dx 



i=Q 



(x 2 — z 2 ) (x — Pi) 

W'(z) (N f + 2) y m (z)P 2{Nc _ r) (z) 



= -( N f + 2)^p + v j : ' m r. zrj;:z<„^,-* ^ 



We used here again the property of contour integration (A. 21). Therefore, we obtain the (A. 20) 
by combining the two contributions (A. 22) and (A. 23) 

TT W '(®cl) = / Tr 1 \ ( w , {z) _ y m (z)P2(N c -r)(z) 



iN c -2N f -A 



z y/P2 (Nc _ r) (z) - 4z 2N f~^A^ 2N f- 4 ' 
Then the second term of above expression can be written as 32 

1 \ y m (z)P2(N c -r)(z) 



Tr- 



Z-Qj ^/P 2 {Nc _ r) (z) - 4 z 2N f -4r+A A AN c -2N f -A 



31 Let us write it as, after an integration over x and using the (A. 11) to change the trace part into 1/pi term 

-§ ^0^ ^^ dx = - LmP ^t Pi) Tr ^ - + 2) L ^ ( a-A a=p<) - For equal massive 

flavors, this becomes — ^ - N ^ + 2^j L ' P2{ ^i~^ ) J^~ Pz " > ■ Due to the mass dependent term containing pi dependent 

term, we can not simplify the second term (A. 20) (as we go to the contour integral around x again, when we 
consider massive flavors) further like as (A. 23) for the massless case. 

32 One can multiply z 2r both in the numerator and denominator and then the orders of the char- 
acteristic polynomial are changed into 2N C . In other words, z 2r P2(N c -r)( z ) = P2N c (z). Remem- 
ber that Tr — h = ^ z -N t +-f™ c ^'} + ^ /+2 ^ and the quantum mechanical expression /Tr— ^ 
± log (P 2{Nc - r) (x) + ^fPi {Nc _ r) (x) - 4x^f-^A^-2N f -^ . 
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Tr- 



y m (z)P 2Nc (z) 



' '(/7$7#) iV/ + 2^ y m (z)P2N c ( z ) 




(7V / + 2 



Therefore, we obtain (A. 20) as 



^W'($ cl )-W'(z) , %T AW'(z)-y m (z)) , . 

''• . , (1 , 7 = + - V . " ^ > */»'< : ; 

Taking into account the relation 

T W'($ d )-W'(z) = /yW'W - W'(z) 



we can write the quantum mechanical expression as follows: 



z-$ \ z - <E» I z - $ 



Tr-^) - ^2)) (W"(s) - y m (*)) (A.24) 

which is the generalized Konishi anomaly equation for the SO(2N c ) case. The resolvent of the 
matrix model R(z) is related to W'(z) —y m (z). When Nf = 0, this result reproduces the result 
in [19] exactly. It does not depend on the r with which we start from. For massive flavors, 
one can generalize the above description similarly, but as we mentioned in the footnote 31, we 
would get rather involved expressions. 

For SO(2N + 1) case, we refer to the version one in the hep-th archive for details. 

Appendix B Addition and multiplication maps 

In [18] it was noted that all the confining vacua of higher rank gauge groups can be con- 
structed from the Coulomb vacua with lower rank gauge groups by using the Chebyshev poly- 
nomial through the multiplication map. Under this map, which was called multiplication map, 
the vacua with the classical gauge group n™=i U (AT,) for a given superpotential are mapped 
to the vacua with the gauge group Hi =1 U(KNi) for the same superpotential where K is a 
multiplication index. This multiplication map was extended to the SO/USp gauge theories 
in [19] and U(N C ) gauge theory with flavors in [20]. In addition to this multiplication map, 
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another map (called the addition map) was introduced in [20]. The addition map reduces many- 
analyses to simpler cases. Thus, we also introduce these two maps, the addition map and the 
multiplication map for the SO(N c ) gauge theory with flavors. 
• Addition map 

Unlike the unitary gauge group, for SO(N c ) with flavors, there are two kinds of addition 
maps: one is for massive flavors and the other is for massless flavors. They should be treated 
separately. As we will show shortly, the addition map with massive flavors can only connect 
the SO(N c ) and SO(M c ) for both N c and M c are an even numbers (or they are odd at the same 
time), but the addition map with massless flavors can connect the SO(2N c ) to SO(2N c + 1). 

Let us start with massive flavors first by using the SO(2N c ) gauge group as an example. 
On assumption that we have two theories, SO(2N c ) with Nf flavors in the r-th branch and 
SO(2N' c ) with Nf flavors in the r'-th branch. On these branches, the SW curves are described 
by 



2 / 2 2\2r 

y 1 = (x — m ) 

„2 _2\2r' 



2 /™\ A„4 \4N'-2N'.-4/ 2 ^2\N'-2r' 



yl = (x a -m a )^[P^_^(x)-4x 4 A^-^/"V-^ ) ' 

For the two curves, if we have the following relations, we find that the polynomials in the 
bracket are identical. 

2(N c -r) = 2(N' c -r'), AN C - 2N f = AN' C - 2N' f , N f - 2r = N' f - 2r'. (B.l) 

These relations are exactly the same as the ones for the U (N c ) case. We expect that certain 
vacua of SO(2N c ) gauge theory with Nf flavors in the r-th branch are related to the one of 
SO(2N' c ) gauge theory with Nf flavors in the r'-th branch. For the SO(2N c + 1) case we have 
the same relation, because the only difference comes from the factor of x 4 A~ 4 which does not 
matter for the relation. They are a common factor in both SW curves. Namely, we have only 
to replace x ±\w c -2N f -i with a .2^4W c -2JV / -2_ It j s a i so noteworthy that starting with SO(N c ) 

for N c even (or odd), we can only get SO(M c ) for M c even (or odd). 

Thus, if on the r-th branch, the SO(N c ) theory with Nf massive flavors has a classical limit 
SO(Nq) x U(Ni) with the effective Nf massless flavors charged under the U(Ni) factor in the 
quartic superpotential, then a SO(N c + 2d) theory with (Nf + 2c?) massive flavors will have a 
classical limit SO(N ) x U(Ni + d) with the effective (Nf + 2d) massless flavors charged under 
U(Ni + d) on the (r + d)-th branch, according to (B.l). Notice that for massive flavors, the 
counting of vacua of the U (Ni) factor is 2N\ — Nf = 2(Ni + d) — (Nf + 2d), which is invariant 
under the above addition map, so it is consistent. In fact, for massive flavors, it is the factor 
U(Ni) with the effective massless flavors so we get the same results as in [20]. 
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Now let us discuss the addition map of massless flavors. For SO(2N c + 1) with N f massless 
flavors, the curve is characterized by 



= p2 Nc{x) _ AK W c -2N f -2 x 2 {x 2 )Nf 

which can be rewritten as 

= P * Nc ( X ) - 4A 4iV c -2(iV / -l)- V(;r 2 ) iV / -l_ (B 2) 

However, the form of (B.2) can be interpreted as the curve of SO(2N c ) with (Nf — 1) massless 
flavors. In other words, SO(2N c + 1) gauge group with Nf massless flavors is equivalent to the 
SO(2N c ) gauge group with (Nf — 1) massless flavors. It is obvious that the above interpretation 
can only be done for massless flavors where SO(2N c ) is related to SO(2N c +1). We call it the 
"special addition maps". In fact, this special addition map comes from a well known result in 
the field theory: giving one flavor in SO(2N + 1) non-zero vacuum expectation value, it can be 
higgsed down to the S0(2N C ) gauge group. 

It seems that except the reinterpretation of S0(2N C + 1) curve to S0(2N C ) curve, we should 
have a similar addition map of massive flavors (we will call it the "general addition map") by 
just setting m = 0. However, we need to be careful in applying the general addition map. As 
we did in the weak and strong coupling analyses, for the SO(N c ) gauge group with massless 
flavors, there is no concept of r-th branch and what we have is the Chebyshev branch or the 
Special branch where the power of x is fixed for a given SO(N c ) factor. The general addition 
map is used to relate the Chebyshev branches (or the Special branches) of two different gauge 
groups. We will give one explicit example to demonstrate the general idea. For SO(2N c ) with 
2M massless flavors, the Chebyshev branch requires a factorization of the curve to be 

V — \ X ) [ r 2N c -(2M+2){ X ) 4A 

Similarly, for SO(2N' c ) with 2M' massless flavors, the curve at the Chebyshev branch is 

y 2 = (* 2 ) 2M ' +2 fe-^ +2) W -4A^- 2 ^'-)] . 



Therefore, if 



N C -(M + 1) = N' C -(M' + 1), 



we can reduce the Chebyshev branch of SO(2N' c ) with 2M' massless flavors to the Chebyshev 
branch of SO(2N c ) with 2M massless flavors. Notice that there is no index r in the above 
relationship. 

• Multiplication map 
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Now we will discuss the multiplication map. The method we used is from the bottom to the 
top, i.e., starting from the known factorization of gauge group with lower rank to the unknown 
factorization of gauge group with higher rank. First let us assume that the factorization of 
curve is given by 



y 



^2N C ( X ) 



2s . AN c -2s-2N 



Ax 2s A 







A(x) = f 2p (x)H 2 N (x 



where s = 2 for the SO(2N c ) group and s = 1 for S0(2N C + 1) group and A(x) the contribution 
of flavors (it can be very general with different masses). Next let us divide out the common 
factor between P 2 2 jv c ( :r ) anc ^ x 2s A(x). For example, if one of the eigenvalues of <3> in the S0(2N C ) 
gauge theory is zero, we can divide out the x A factor at the two sides of the curve. Or if A(x) 
has factor (x 2 — m 2 ) 2 and we are at the r = 1 branch, the factor (x 2 —m 2 ) 2 can also be divided 
out. After dividing out all the common factors, the remaining curve is given by 



P 2 M (x) - 4At Nc ~ 2s - 2Nf A(x) = f^(x)H 2 M _~(x) 



where p is some number. Now let us define 



Pm(x) 



X 



2 V A 



2N c -s-N f 



(B.3) 



and 



Pk uU) = //Af ^ X \~H^Y T K (r) 



(B.4) 



with rf K = 1. It is worth to notice that although x has \J A(x) in the denominator, the Prm{x) 
is a perfect polynomial of x. The reason is very simple: Txit) has only an even (odd) power 

K—q „ Kq 

of t in the polynomial if K is an even (odd) number, so y A{x) = A(x)^~ is a polynomial 
of x because (K — q) is always an even number 33 . Using (B.3) and (B.4) it is easy to see that 
one obtains 



P 2 KM {x)-Ak^- 2s ' 2N ^ K A{x) K 

( 

^{W c -2s-2N f )K A( X ) K 



Pkm(x) 



x 



= A^" 2s ~ 2N ^ K A{x) K \T 2 {x) - 1 



2 ( ^-^^Mx^Uk^x) (x 2 - l) 



K 



— 1 



33 We want to thank Freddy Cachazo and Oleg Lunin for discussing this general situation of the multiplication 
map. 
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K-l 



Notice that a similar reason guarantees that y A(x) k ~ 1 Uk-i(x) is a polynomial of x. The above 
calculation shows that we have the solution of factorization of Pkm(x) with matter dependent 
part A(x) K . Finally to recover to the SO(M c ) gauge group, we need to multiply back some 
factor of x p to make sure that (KM + p/2) is even and x p A(x) K has x 2 or x 4 factor. 
To demonstrate the above general results, let us consider several concrete examples: 

1) . The first example is that we consider the breaking pattern SO(2N c ) — > ]T?=i U(Nj) at 
the r = branch, so the P2N C ( X ) does not have a common factor with x 4 A(x) = x A (x 2 — m 2 ) N f. 
34 Putting M = 2N C and A(x) = x 4 (x 2 -m 2 ) N f , we get A(x) K = (x 2 -m 2 ) KN fx 4K . This means 
that the multiplication map lifts the original solutions to SO(2KN c ) with KNf massive flavors 
and 2(K — 1) massless flavors (or SO(2KN c + 1) with KNf massive flavors and 2(K — 1) + 1 
massless flavors). It is somewhat surprising that we get the massless flavors after lifting. The 
reason is that the original breaking SO(2N c ) — > Ylj =1 U(Nj) can be effectively considered as 
U(N C ) — > YYj = iU(Nj) with Nf massive flavors plus two massless flavors which come from a 
factor x 4 = (x 2 ) 2 by looking at the SW curve. From this new point of view, we have massless 
flavors at the beginning. 

2) . For the breaking SO(2N c ) -> SO(2N ) x U.]=iU(Nj) : we can divide out at least the 

factor x A from the SW curve because the characteristic function P2N C { X ) an d 4x 4 Aq JVc 4 2Nf A(x) 

possess a factor x 2 , so M = (2N c — 2) and A(x) = (x 2 — m 2 ) N f . The multiplication map will give 

the factorization of SO(K(2N c — 2) + 2) with KNf flavors. It is significant that the +2 of the 

rank K(2N c -2) +2 comes from putting back the x 4 factor to P KM (x) -AK^ Nc ~ 2s - 2N ^ K A(x) K . 

Since we can not multiply a factor x 2 back ((KM + p/2) is not even), it is impossible to obtain 

the map from SO(2N c ) to SO(2N c + 1) which reveals the same phenomena observed in [19] 

for the pure SO(2N c ) gauge theory. However, if there are massless flavors, one gets from 

SO(2N c ) to SO(2M c + 1) through the special addition map. The multiplication map has 

another application. We have seen that the relation between f2 P (x) and W'(x) is changed in 

certain situation under the presence of the flavors. For SO(2N c ) gauge theory with Nf flavors 

at the n-th branch, when n > —Nf + 2N c — 3, the relation between f2 P (x) and W'(x) is modified. 

On the other hand, for SO(2KN c — 2K + 2) gauge theory with KNf flavors at n-th branch, 

this happens when n > —KNf + (2KN C — 2K + 2) — 3. Since we are considering the particular 

case Nf < 2N C — 2, for asymptotic free theory, the relation between f2 P (x) and W'(x) will be 

34 For simplicity, we have assumed that all Nf flavors have the same mass. Generalization to different masses 
is obvious. 
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different for SO(2N c ) gauge and SO(2KN c — 2K + 2) gauge theories. If K is large enough, the 
presence of flavors will not modify the relationship between f2 P {x) and W'(x) and the geometric 
proofs [10, 11, 12] can go through [20]. In other words, the geometric picture is really for a 
large N c limit. 

3). To break SO(2N c +l) -> SO(2N + 1) xf]J =1 U(Nj), we can divide out at most the factor 
x 2 if all masses are not zero, so M = (2N C — 1) and A(x) = (x 2 — m 2 ) N f . The multiplication map 
will give us KM = K(2N C — 1). Notice that when we multiply back x p we require (KM +p/2) 
is even. If K is even, we can only multiply back x 4 factor so it is SO(K(2N c — 1) + 2) with 
KNf flavors. If K is odd, we can only multiply back x 2 factor so it is SO(K(2N c — 1) + 1) 
with KNf flavors. Notice that this conclusion is made under the assumption that there are no 
massless flavors at the beginning. If there are massless flavors, by the "special addition maps" , 
we can jump from SO(2N c ) gauge theory to SO(2N c + 1) gauge theory. 
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